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ABSTRACT 


Pending  changes  in  NRC  emergency  planning  requirements  have 
prompted  this  reevaluation  of  the  maximum  radiation  dose  to  an 
individual  located  at  the  exclusion  area  boundary  of  the  MIT  Reactor 
during  the  first  two  hours  of  its  Design  Basis  Accident.  The  DBA  is 
postulated  to  be  the  melting  of  four  fuel  plates  in  one  element, 
releasing  a  maximum  of  1.5  %  of  the  core  fission  product  inventory 
to  the  coolant. 

The  approach  used  was  to  evaluate  the  major  release  pathways  to 
the  exclusion  area  boundary.  First  the  reactor  system  was  analyzed  to 
determine  the  fission  product  release  fractions  to  the  containment. 

The  dose  due  to  leakage  from  the  containment  was  evaluated  using  a 
standard  Guassian  diffusion  model  with  local  meteorological  data. 

Gamma  radiation  reaching  the  boundary  by  simple  direct  penetration  of 
the  containment  was  determined  from  standard  shielding  manual 
approaches.  A  compton  scattering  model  was  developed  and  applied  tc 
photons  scattering  from  air  (skyshine)  and  from  the  steel  containment 
roof.  Finally,  the  largest  containment  penetration,  a  truck  airlock, 
was  checked  for  radiation  streaming. 

The  resulting  total  maximum  dose  is  estimated  to  be  0.595  rads  to 
the  whole-body  and  0.118  rads  to  the  thyroid.  These  are  well  within 
the  regulatory  limits  for  exclusion  area  releases  during  an  accident. 
The  steel  scattering  dose  was  the  largest  component  of  the  total  whole- 
body  dose  (63%)  and  the  truck  lock  penetration  dose  was  found  to  be 
negligible.  In  addition,  the  behavior  of  each  component  of  the  dose 
beyond  the  exclusion  area  boundary  was  estimated  based  on  a  few 
prominent  gamma  energies  with  the  resulting  conclusion  that  the  total 
dose  decreases  beyond  the  boundary. 
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ABSTRACT 

Pending  changes  in  NRC  emergency  planning  requirements  have 
prompted  this  reevaluation  of  the  maximum  radiation  dose  to  an 
individual  located  at  the  exclusion  area  boundary  of  the  MIT  Reactor 
during  the  first  two  hours  of  its  Design  Basis  Accident.  The  DBA  is 
postulated  to  be  the  melting  of  four  fuel  plates  in  one  element, 
releasing  a  maximum  of  1.5  %  of  the  core  fission  product  inventory 
to  the  coolant. 

'The  approach  used  was  to  evaluate  the  major  release  pathways  to 
the  exclusion  area  boundary.  First  the  reactor  system  was  analyzed  to 
determine  the  fission  product  release  fractions  to  the  containment. 

The  dose  due  to  leakage  from  the  containment  was  evaluated  using  a 
standard  Guassian  diffusion  model  with  local  meteorological  data. 

Gamma  radiation  reaching  the  boundary  by  simple  direct  penetration  of 
the  containment  was  determined  from  standard  shielding  manual 
approaches.  A  compton  scattering  model  was  developed  and  applied  to 
photons  scattering  from  air  (skyshine)  and  from  the  steel  containment 
roof.  Finally,  the  largest  containment  penetration,  a  truck  airlock, 
was  checked  for  radiation  streaming. 

The  resulting  total  maximum  dose  is  estimated  to  be  0.595  rads  to 
the  whole-body  and  0.118  rads  to  the  thyroid.  These  are  well  within 
the  regulatory  limits  for  exclusion  area  releases  during  an  accident. 
The  steel  scattering  dose  was  the  largest  component  of  the  total  whole- 
body  dose  (63%)  and  the  truck  lock  penetration  dose  was  found  to  be 
negligible.  In  addition,  the  behavior  of  each  component  of  the  dose 
beyond  the  exclusion  area  boundary  was  estimated  based  on  a  few 
prominent  gamma  energies  with  the  resulting  conclusion  that  the  total 
dose  decreases  beyond  the  boundary, 
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CHAPTER  I 


INTRODUCTION 

1.1  Purpose 

The  purpose  of  this  work  is  to  determine  a  conservative  estimate 
of  the  maximum  amount  of  radiation  an  individual  located  at  the 
exclusion  area  boundary  of  the  MIT  research  reactor  (MITR-II)  would 
receive  during  the  first  two  hours  of  the  reactor's  design  basis 
accident  (DBA) .  The  motivation  for  this  work  stems  from  new  guidelines 
forthcoming  from  the  Nuclear  Regulatory  Commission  (NRC)  requiring 
the  establishment  of  comprehensive  emergency  plans  for  accident 
response  by  each  research  reactor  facility.  The  classification  of 
the  emergency  and  the  appropriate  actions  are  related  to  the  size 
of  the  radiation  release. 

The  MIT  reactor  is  of  particular  interest  as  it  is  located  in 
the  City  of  Cambridge  within  the  limits  of  metropolitan  Boston. 

The  facility  consists  of  a  containment  structure,  cooling  towers 
and  a  support  building  and  is  located  in  an  industrial  area  just 
near  the  main  MIT  campus.  The  surrounding  buildings  are  normally 
occupied  during  the  work  week  and  the  street  directly  in  front  of 
the  reactor  contains  a  heavy  flow  of  vehicular  and  pedestrian  traffic. 

The  exclusion  area  consists  of  the  reactor  parking  lot,  cooling  tower 
area  and  a  portion  of  the  adjacent  building  NW-12  as  shown  in  Fig.  1.1[1] 

1.2  Description  of  the  MITR-II 

The  present  MIT  reactor  (MITR-II)  is  a  heavy-water  reflected,  light- 


Fipure  ].]  Exclusion  Area  Boundaries  [1] 
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water  cooled  and  moderated  5  megawatt  research  reactor.  It  is 

fueled  by  93%  enriched  U-Al^  flat-plate  fuel  elements  arranged  in 

a  compact  hexagonal  core.  This  core  design  maximizes  the  neutron 

flux  in  the  D  0  reflector  region  where  numerous  experimental  beam 
2 

ports  are  located.  The  core  is  contained  within  a  light-  water  filled 
aluminum  tank  which  is  in  turn  contained  within  the  D2O  reflector 
tank.  The  H^O  coolant  is  directed  so  as  to  flow  down  along  the  tank 
walls  and  then  upwards  through  the  fuel  elements.  Heat  from  the 
primary  systems  is  transferred  by  heat  exchangers  to  the  secondary 
system  which  dissipates  it  to  the  atmosphere  through  the  cooling 
towers . 

The  reactor  is  located  at  the  center  of  a  gastight  cylindrical 
steel  building  equipped  with  a  controlled  pressure  relief  system. 

Access  to  the  reactor  building  is  through  either  of  two  personnel  or  one 
truck  air-locks.  Al)  building  penetrations  are  either  sealed  perm¬ 
anently  or  can  be  sealed  rapidly  by  manual  or  automatic  operation. 

The  building  is  designed  to  withstand  a  maximum  overpressure  of 
2  psi  and  normally  operates  at  a  slight  negative  pressure.  The  shell  has 
0.95  cm  (3/8  in)  thick  steel  plates  on  the  sides  (5/8  in  on  the  dome) 
with  an  outside  diameter  of  22.6  ra  (74  ft)  and  a  height  of  14.9  m  (49  ft). 
Contained  within  the  circumference  of  the  steel  shell  is  a  cylindrical 
concrete  shadow  shield  0.61  m  (2  ft)  thick  and  9.6  m  (31.5  ft)  high. 

For  a  more  complete  description  of  the  reactor  and  facility  the  reader 
is  referred  to  the  Reactor  Systems  Manual  [1], 
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1.3  The  MITR-II  Design  Basis  Accident 

The  design  basis  accident  is  the  maximum  credible  accident 
which  could  result  in  the  release  of  radiation  from  the  reactor. 

For  the  MITR-II,  this  is  postulated  to  be  a  coolant  flow  blockage 
in  the  hottest  channel  of  the  center  fuel  element.  This  could  occur 
as  the  result  of  some  foreign  material  falling  into  the  reactor  during 
refueling.  After  the  pumps  are  started,  the  material  would  be  swept 
from  the  bottom  of  the  tank  up  to  the  fuel  element  flow  entrance. 

Due  to  the  size  of  the  openings  in  the  adapters  at  the  end  of  each 
fuel  element  any  piece  of  material  which  could  pass  through  the 
adapters  would  not  be  large  enough  to  restrict  the  flow  in  more 
than  five  coolant  channels  of  one  element,  leading  to  overheating 
in  a  maximum  of  four  fuel  plates.  It  is  conservatively  assumed  that 
the  entire  active  portions  of  all  four  plates  melt  completely,  releasing 
their  inventory  of  fission  products  to  the  coolant  water. 

1.4  Organization  of  this  Thesis 

Since  the  MITR-II  is  not  a  closed  system,  it  is  possible  for 
fission  products  in  the  coolant  to  emerge  directly  into  the  contain¬ 
ment  atmosphere.  The  quantity  of  fission  products  in  the  melted 
fuel  and  the  fraction  which  is  released  to,  and  remains  airborne 
in,  the  containment  will  be  determined  in  Chapter  2.  The  airborne 
fission  products  are  subject  to  release  from  the  containment  by  leakage 


through  small  cracks  and  gaps  and  by  venting  through  the  pressure 
relief  system.  The  dose  expected  from  these  releases  is  investigated 
in  Chapter  3.  Those  fission  products  which  do  not  escape  the  contain¬ 
ment  can  still  add  to  the  dose  at  the  exclusion  area  boundary  as  a 
source  of  gamma  radiation.  This  can  reach  the  boundary  as  simple 
unattenuated  direct  radiation  through  the  containment  sides  or  as 
scattered  radiation  which  reaches  the  ground  through  interactions 
with  the  air  or  the  steel  containment  roof.  The  direct  dose  is 
treated  in  Chapter  4  and  the  scattered  dose  in  Chapter  5.  Finally, 
the  effect  of  gamma  release  through  containment  penetrations  is 
investigated  in  Chapter  6  by  examination  of  the  truck  lock. 

1.5  Previous  Work 

The  first  estimate  of  an  exclusion  area  radiation  release  was 
made  in  the  Final  Hazards  Summary  Report  [2]  for  the  original 
MIT  reactor,  MITR-I.  This  was  a  five  megawatt  heavy-water  moderated 
and  cooled  reactor  located  in  the  same  building  and  site  as  the 
present  MITR-II.  Although  the  DBA  and  fission  product  release  for  the 
MITR-I  were  different,  the  report  estimated  the  exposure  due  to 
leakage,  direct  and  scattered  radiation  and  formed  the  basis  for 
all  subsequent  analyses. 

Prior  to  the  reactor  being  converted  to  light-water  coolant 
(along  with  other  changes)  in  1974  a  Safety  Analysis  Report  (SAR)  [3] 


was  published  which  proposed  the  present  DBA  and  updated  the  estimate 
of  exclusion  area  release.  The  fission  product  inventory  and  release 
estimates  were  based  heavily  on  TID-14844  [4]  while  the  dose  calcul¬ 
ations  were  very  similar  to  the  Final  Hazards  Report.  The  SAR  intro¬ 
duced  an  added  dose  due  to  the  pressure  relief  system  and  modified 
the  leakage  and  steel  scattering  calculations. 

The  latest  work  is  a  thesis  by  McCauley  [5].  In  it,  he  makes 
a  convincing-  argument  that  the  present  DBA  is  in  fact  the  maximum 
credible  accident.  In  addition,  he  revises  the  release  fractions, 
particularly  those  of  iodine,  based  on  more  recent  studies  and 
experience  (such  as  the  Three  Mile  Island  accident) . 

Although  the  purpose  of  this  thesis  is  to  calculate  the  same 
doses  via  the  same  release  paths, the  approach  will  differ  significantly 
The  fission  product  source  term  used  in  the  SAR  includes  only  three 
elements;  Xenon,  Krypton,  and  Iodine.  In  Chapter  2,  the  entire 
spectrum  of  possible  fission  products  will  be  considered.  The 
calculations  of  the  leakage  dose  is  revised  in  accordance  with 
recently  published  NRC  atmospheric  diffusion  model  requirements. 

The  direct  dose  component  more  accurately  takes  into  account  the 
shape  and  size  of  the  volumetric  source.  By  far,  the  weakest  portion 
of  the  SAR  is  that  concerning  air  scattering.  Due  to  its  grossly 
conservative  approach,  the  SAR  predicts  air  scattering  to  be  the 
dominant  component  of  the  total  dose.  In  Chapter  5  a  more  exact 
scattering  model  is  developed,  which,  while  still  conservative,  is 
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CHAPTER  2 

DEVELOPMENT  OF  THE  CONTAINMENT  SOURCE  TERM 


2.1  Fission  Product  Build-up  in  the  Fuel 

The  fission  product  inventory  in  the  fuel  at  the  time  of  the 

accident  is  assumed  to  be  equal  to  the  maximum  value  of  equilibrium 

fission  products  after  an  essentially  indefinite  period  of  irradiation 

at  5  MW.  This  is  consistent  with  the  present  SAR[3]and  NRC  guidelines 

for  exclusion  area  calculations  as  presented  in  TID-14844 [4] .This  is 

a  conservative  assumption  in  that  the  present  five  days  per  week  operating 

schedule  does  not  result  in  saturation  of  most  important  fission 
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products,  including  I  [5].  However,  if  the  operating  schedule 
were  to  be  expanded  at  some  future  time,  the  assumption  of  equilibrium 
fission  product  inventories  could  become  more  realistic. 

The  fission  product  isotopes  of  interest  have  been  distilled  from 
a  suggested  list  in  Thompson  and  Beckerley[6]and  from  those  used 
in  the  Reactor  Safety  Study,  WASH-1400[7] .  The  selection  has  been  based 
on  volatility,  quantity  produced,  half-life  and  degree  of  biological 
effectiveness . 

Highly  volatile  isotopes  are  of  primary  concern  as  they  can 
readily  escape  to  the  containment  atmosphere.  Low  volatility 
isotopes  (the  "solids")  are  of  less  importance  but  still  have  some 
finite  escape  probability.  Isotopes  which  are  not  produced  in  sufficient 
quantity  to  have  a  measurable  release  can  be  neglected.  Isotopes  with 
half-lives  of  less  than  ten  minutes  have  been  neglected  since  they  would 
be  effectively  gone  within  the  first  half  of  the  period  of  the  accident. 
Conversely,  some  isotopes  which  were  not  considered  in  WASH-1400  due 
to  their  short  half-lives  (with  respect  to  the  time  scale  of  that  study) 
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have'  been  included  in  this  study.  Finally,  only  isotopes  which 
contribute  significantly  to  either  whole-body  or  thyroid  doses  have 
been  considered.  The  criteria  for  inclusion  was  a  biological  effective 
ness  within  three  orders  of  magnitude  of  the  most  effective.  The 
resulting  isotopes  are  listed  in  Table  2.1. 

The  saturation  activities  of  the  fission  products  have  been 
determined  using  two  methods,  analytical  and  computational .  Each 
method  is  discussed  in  the  following  sections. 


E 


so  that. 


V -Z,(J)  Y. 
f  f  I  1 


where  N^  =  Saturation  number  of  nuclei  of  isotope  i. 

This  can  be  expressed  as  a  function  of  reactor  power  using  the  fact 
that  one  megawatt  is  produced  by  3.2  x  1016  fissions/sec,  therefore 


Vfr^  =  P ( 3 . 2  x  10  )  fissions/sec. 

where,  P  =  Reactor  power  (MW), 

„,i  P(3.2  x  1016)  Y-r 

and  consequently  N  =  - ^ - - - — 

S  A  . 


(2.2) 


By  definition,  the  activity  (in  disintegrations/sec)  of  is 

Since  one  curie  equals  3.7  x  10^  dis/sec,  the  saturation  activity, 

Q^,  due  to  the  presence  of  N*,  is 

N1X, 

O1  _■  s  1 

^s  3.7  x  10-^  curies 


,i  P  (3 . 2  x  10iD)  Y,- 
}s  ~  3.7  x  10i0 - 


resulting  in 


Q1  =  8.65  x  10  5  Py.  curies 
s  1 


(2.3) 


Note  that  the  above  expression  implicitly  assumes  a  constant  rate  of 
production  of  isotope  i  (i.e.  no  fuel  depletion).  It  also  does  not 
take  into  account  the  loss  or  addition  of  fission  products  due  to 
neutron  capture.  Since  most  fission  products  are  formed  by  decay 


of  radioactive  parents  as  well  as  direct  production  by  fission,  the 
yield  value  used  must  be  the  total  yield  due  to  both  effects.  Values 
of  Y4  and  for  each  isotope  are  listed  in  Table  2.1. 
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TABLE  2.1 

Total  Core  Fission  Product  Inventory 
(in  order  of  decreasing  volality)  [6,7,8] 


Isotope 

Half-life 

Decay 

Total 

Saturation 

Cl/2 

Cons  tant 

Xh 

(sec  X) 

Yield 

Yi 

(%) 

Activity 

«s  + 
(Ci) 

Kr 

85m 

4. 36h 

4.41  x  10"5 

1.5 

6.49  x  104 

87 

78m 

1.48  x  10"4 

2.7 

1.17  x  105 

88 

2 . 77h 

6.95  x  10"5 

3.7 

1.60  x  105 

Xe 

131m 

12.  Od 

6.68  x  10“ 7 

0.03 

1.30  x  103 

133m 

2.3d 

3.49  x  10“6 

0.16 

6.92  x  103 

133 

5 . 27d 

1.52  x  10“6 

6.5 

2.81  x  105 

135m 

15.6m 

7.40  x  10~4 

1.8 

7.78  x  104 

135 

9 . 13h 

2.11  x  10“5 

6.2 

4.13  x  104 

138 

17m 

6.79  x  10“4 

5.5 

2.38  x  105 

I 

131 

8.05d 

9.96  x  10“ 7 

2.9 

1.25  x  105 

132 

2 . 4h 

8.02  x  10“5 

4.4 

1.90  x  105 

133 

20. 8h 

9.25  x  10“6 

6.5 

2.81  x  105 

134 

52.5m 

2.20  x  10“5 

7.6 

3.29  x  105 

135 

6.68h 

2.89  x  10“5 

5.9 

2.55  x  105 

Br 

83 

2 . 4h 

8.02  x  10“5 

0.48 

2.08  x  104 

84 

30m 

3.85  x  10“4 

1.1 

4.76  x  104 

Cs 

134 

2.0v 

1.10  x  10" 8 

•k 

0.0 

2.86  x  105 

136 

13d 

6.17  x  10" 7 

0.006* 

4.14  x  104 

137 

26. 6y 

8.27  x  10“10 

5.9 

2.31  x  105 

Rb 

86 

19. 5d 

4.11  x  10" 7 

-5* 

2.8x10 

6.12  x  104 

x  -  Analytical  value  unless  otherwise  noted 
t  -  Computational  value  *  -  Direct  yield 


S-t  -H-  -H-  -H- 
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TABLE  2.1  (Cont.) 


Isotope 

Half-lif e 

Decay 

Constant 

Total 

Yield 

Te  127m 

90d 

8.82  x  10"8 

0.056 

127 

9 . 3h 

2.07  x  10“5 

0.25 

129m 

33d 

2.43  x  10"7 

0.34 

129 

72m 

1.60  x  1 o“4 

1.0 

131m 

30h 

6.42  x  10-3 

0.44 

131 

24 . 8m 

4.66  x  10~4 

2.9 

132 

77h 

2.50  x  10~6 

4.4 

133m- 

63m 

1.83  x  10~4 

4.6 

134 

44m 

2.63  x  10“4 

6.7 

Sr  91 

9  7h 

1.99  x  10"5 

5.9 

Ba  140 

12. 8d 

6.27  x  10-7 

6.3 

Ru  103 

41d 

1.96  x  10-7 

2.9 

105 

4 . 5h 

4.28  x  10“3 

0.9 

10  6 

l.Oy 

2.20  x  10~8 

0.38 

Rh  105 

36. 5h 

5.27  x  10"6 

0.9 

Tc  99m 

6.04h 

3.19  x  10~3 

~0.6 

Mo  99 

6  7h 

2.88  x  10~6 

6.1 

Sb  127 

9  3h 

2.07  x  10'6 

0.25 

129 

4.6h 

4.32  x  10~5 

1.0 

Nd  147 

11.3d 

7.10  x  10“ 7 

2.6 

La  140 

40. 2h 

4.79  x  10“6 

6.3 

Ce  141 

32d 

2.51  x  10“7 

6.0 

143 

32h 

6.01  x  10-6 

6.2 

144 

290d 

2.76  x  10“8 

6.1 

Zr  95 

63d 

1.27  x  10"7 

6.4 

97 

17h 

1.13  x  10-3 

6.2 

Nb  95 

35d 

2.29  x  10“7 

6.4 

Saturation 
Activi ty 

2.42  x  103 
1.08  x  10  4 
1.47  x  104 
<4.32  x  104 

1.90  x  104 
1.25  x  105 

1.90  x  105 
1.99  x  105 

2.90  x  105 

2.55  x  105 

2.72  x  105 

1.25  x  105 
3.89  x  1 04 

1.64  x  104 

3.89  x  104 
2.59  x  104 

2.64  x  105 

1.08  x  104 
4 . 32  x  105 

1.12  x  105 

2.72  x  105 

2.59  x  105 
2.68  x  105 

2.64  x  105 

2.77  x  105 
2.68  x  105 

2.77  x  105 
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2.1. '2  Computational  Calculation 


The  analytical  method  cannot  be  used  to  determine  the  (Rvalues  of 
all  isotopes  since  some  isotopes  have  significant  amounts  of  production 
and  reduction  due  to  neutron  capture.  Values  of  for  these  isotopes 
(Xe  135,  Cs  134,  136,137,  and  Rb  86)  can  be  estimated  using  computational 
results. 

A  computer  code  to  calculate  fission  product  production  has  been 
developed  by  Blomeke  and  Todd[8]which  solves  eleven  simultaneous 
equations  of  fission  product  build-up  and  decay.  This  model  also 
assumes  no  depletion  of  the  fuel  but  includes  neutron  capture  reactions. 
Values  of  N.  for  various  combinations  of  thermal  flux,  irradiation  time 

l 

and  decay  time  have  been  calculated  and  are  available  in  graphic 

form.  For  this  study  only  the  saturation  values  are  required.  These 

are  expressed  as  the  saturated  number  of  fission  product  atoms 

235  i  ° 

produced  per  initial  atom  of  U  (N  /N„,-)  as  a  function  of  thermal 

S  « 

13  Z 

flux.  The  core  averaged  thermal  flux  value  of 'V 4  x  10  neutrons/cm 

i  ° 

sec  has  been  used  to  enter  the  graphs.  The  resulting  values  of  N  /N„r 

s  25 

for  the  neutron-capture  influenced  isotopes  are  displayed  in  Appendix  A, 
Table  A.l. 


The  resulting  activity  of  each  isotope  can  be  determined  without 


actual  knowledge  of  .  As  before,  the  saturation  activity  is 

,i 


i 

^s  = 


A.N 

1  s 


3.7  x  10 


10 


curies 


or 


«> 


<X.)(n”5>  <">25 > 


3.7  x  10 


10 


curies 


25 


>16, 


25 


and  given  P(3.2  x  10  )  =  Vf£f4>T  =  N25<jf  <J>T 


„°  .  ,  „°  PC3.2  x  1016) 

can  be  expressed  as  - —  atoms 

25 

where  =  Microscopic  fission  cross-section 

O  c  t 

for  U  =  580  x  10“  cm 


(2.5) 


Substituting  for  N  the  saturation  activity  becomes 


1.49  x  1027  X.P(N1/N°  ) 

Q1  =  - -r-i - 

s  <|>_ 


curies 


(2.6) 


25  “ 2  A  2 

(Note:  The  2200  m/s  value  0^  =  580  x  10  cm  without  temperature 

adjustment  must  be  used  to  be  consistent  with  the  Blomeke  and  Todd 
derivation) . 

Values  of  for  the  above  named  istotopes  obtained  using  Eq .  2.6  are 
listed  in  Table  2.1. 

As  a  check  of  the  analytical  method,  values  of  have  been 
calculated  using  Eq.  2.6  for  all  the  non-neutron-capture  influenced 
isotopes.  In  all  cases  the  analytical  and  computational  values 
agree  within  slight  variances  due  to  graphical  interpolation.  Because 
of  this  need  to  interpolate  in  the  Blomeke  and  Todd  method,  the 
analytical  values  are  considered  more  accurate  and  have  been  used 
in  all  subsequent  calculations  for  the  non-neutron-capture  influenced 
isotopes . 


2.1.3.  Fission  Product  Inventory  in  the  Melted  Fuel 


In  the  previous  sections  the  saturated  core  inventory  of  fission 
product  activities  has  been  determined.  Only  a  small  portion  of  this 
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core'  inventory  will  be  contained  in  the  four  fuel  plates  which  are 
assumed  to  melt.  If  the  core  contains  25  fuel  elements  and  the  radial 
peaking  factor  is  1.45  [3]  then  the  fraction  of  the  total  saturated 
core  inventory  which  is  contained  in  the  four  center  fuel  plates,  F  , 


can  be  determined  to  be 


s 


F  =  _ (4  plates)  _  (1,45) _  =  Q15 

s  (15  plates/element)  (25  elements) 

Therefore,  a  maximum  of  1.5%  of  each  is  available  for  release 
from  the  melted  fuel. 


2.2  Release  Fractions 

2.2.1  Overview  of  Release  Fractions 

While  the  estimation  of  fission  product  accumulation  can  be 
done  with  reasonable  accuracy,  that  of  release  fractions  cannot. 
Although  numerous  experiemnts  and  calculations  have  been  performed 
for  many  combinations  of  temperature,  atmosphere,  fuel  types  and  the 
like,  there  is  no  universal  agreement  as  to  what  the  exact  release 
fraction  for  any  particular  accident  will  be.  At  best,  only  an 
estimate  of  the  range  of  possible  release  fractions  can  be  made. 

For  estimates  of  potential  releases  of  radioactivity  to  the  public 
a  conservative  bias  is  desired.  Although  not  strictly  binding, 

10  CFR  100,  "Reactor  Site  Criteria"! 9 ]  suggests  use  of  the  release 
fractions  in  TID-14844,  which  are:  100%  of  the  noble  gases,  50%  of 
the  halogens  and  1%  of  the  solids  released  to  the  containment;  and 
of  the  iodine  released  to  the  containment,  50%  remains  available  for 
release  from  the  containment  [4]. 
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NRC  Regulatory  Guides  1.3  and  1.4  require  the  TID  source  term 
to  be  used  in  evaluating  the  consequences  of  BWR  and  PWR  LOCAs,  with 
the  exception  that  the  release  of  solids  is  neglected  [10,11]. 

Tha  above  are  general  estiamtes  not  particularly  suited  to  the 
MIT  reactor.  In  the  following  sections  a  conservative  estimate 
more  applicable  to  the  MITR  will  be  made. 

2.2.2  Release  from  the  Fuel 

Studies  conducted  on  the  melting  of  uranium  metal  and  uranium- 

aluminum  show  that  the  quantity  of  fission  products  released  increases 

with  increasing  burn-up,  percentage  of  melting,  temperature  of  melting 

and  the  total  time  of  melting  1.6,12,13,14] .  The- use  of  the  data  is 

limited  by  the  fact  that  it  is  given  for  atmospheres  of  air,  steam- 

air  and  helium  only.  Assuming  that  the  melting  of  the  fuel  is 

accompanied  by  film  boiling,  data  for  steam-air  mixtures  can  be  used. 

Assuming  also  that  the  fuel  plates  have  40%  burn-up  and  are  completely 

melted  at  a  temperature  of  1100°C  (2012°F)  over  a  period  of  60  minutes, 

the  percentage  of  fission  products  contained  in  the  four  melted  plates 

which  are  released  from  the  fuel,  ,  can  be  conservatively  estimated  as 

100%  of  the  noble  gases  (Kr,Xe) 

100%  of  the  halogens  (I,  Br) 

70%  of  the  Tellurium 

30%  of  the  alkali  metals  (Cs ,  Rb) 

1%  of  the  remaining  fission  products 

These  values  are  comparable  to  the  meldown  release  component  values 
used  in  WASH-1400  [  15  ]  .The  differences  can  be  attributed  to  the  difference 
in  fuel  (UO^  clad  with  zirconium  vs.  U-Al)  and  the  type  of  accident 


(core  meltdown  with  loss  of  coolant  vs.  fuel  plate  melting  with  primary 
envelope  intact). 

2.2.3  Release  From  the  Primary  System 

The  two  methods  of  fission  product  depletion  in  the  primary 
system  are  deposition  on  core  surfaces  (plateout)  and  absorption  in 
the  coolant.  Deposition  is  the  dominant  mechanism  for  loss-of-coolant 
accidents,  and  retention  in  the  coolant  applies  to  accidents  in  which 
the  core  remains  covered.  Most  research  has  involved  the  first  case, 
since  these  accidents  contribute  the  greatest  amount  of  public  risk. 

The  limited  work  in  the  second  area,  as  reported  in  WASH-1400 [ 15 ] . 
indicates  a  100%  release  of  the  noble  gases  through  the  coolant,  and 
a  range  of  100%  to  0.1%  for  the  remaining  isotopes. 

The  retention  of  iodine  has  received  considerable  attention  in 
the  last  two  years.  Recent  studies  indicate  that  iodine  may  be  more  like¬ 
ly  to  be  released  from  the  fuel  in  the  form  of  Csl  rather  than 
elemental  iodine [6].  The  importance  of  this  is  that  Csl  is  extremely 
soluble  in  water  and  can  be  expected  to  be  retained  to  a  high  degree. 

A  computer  simulation[16]has  shown  that  the  release  fraction  for 
through  water  is  33%  while  that  of  Csl  is  1%.  The  formation  of  Csl  is 
enhanced  by  a  reducing  atmosphere,  w’hile  in  the  presence  of  oxygen 
(including  steam,  but  particularly  air)  Csl  rapidly  breaks  down  resulting 
in  the  formation  of  I  or  I  .  Although  one  can  argue  over  which  form 
of  iodine  will  be  produced  by  UO^ ,  it  appears  clear  that  low-oxygen  U-Al 
fuel  will  produce  Csl.  McCauley[5]shows  that  sufficient  Cs  will  exist 
to  combine  with  all  available  iodine  and  that  all  the  Csl  will  be 
dissolved  in  the  coolant.  In  addition,  actual  measurements  of  the 
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activity  dua  to  I  in  the  coolant  and  that  released  to  the  exhaust 
filters  during  normal  operation  indicates  a  release  fraction  of  0.0l%[5] 
Retention  of  the  remaining  fission  products  in  water  has  not  been 
investigated  in  the  literature  but  can  be  expected  to  be  low.  Reten¬ 
tion  of  all  fission  products  will  be  primarily  a  function  of  their 
solubility  and  vapor  pressure,  which  will  change  depending  on  the 
thermodynamic  conditions  in  the  reactor  vessel.  The  above  dicussion 
notwithstanding,  conservative  estimates  of  the  fraction  of  fission 
products  released  from  the  fuel  which  escape  the  primary  system,  F  , 

p 

will  be  assumed  to  be: 

100%  of  the  noble  gases 
10%  of  all  other  isotopes. 


2.2.4  Natural  Depletion  in  Containment 

As  the  MITR  has  no  containment  sprays  or  other  engineered  safety 

features  to  reduce  the  quantity  of  fission  products  in  the  containment 

atmosphere,  depletion  of  the  radioactive  isotopes  released  to  the 

containment  can  occur  only  through  natural  processes.  These  include 

diffusion  to  and  deposition  on  the  walls  and  equipment,  agglomeration 

and  gravitational  deposition  on  the  floor,  and  condensation  in  the 

steam.  The  fraction  of  fission  products  released  to  the  containment 

which  remain  airborne  in  the  containment  atmosphere  will  be  designated 

as  F  . 
c 

Tne  noble  gases  are  not  expected  to  be  depleted  by  any  of  the 
above  methods  resulting  in  a  release  fraction  of  100%. 
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Iodine  will  be  assumed  to  be  present  in  elemental  form.  Even 
if  released  as  Csl,  the  highly  oxidizing  atmosphere  in  the  containment 
will  cause  it  to  rapidly  dissociate.  For  iodine  the  primary  method 
of  depletion  is  deposition  on  the  walls  through  natural  convection. 


For  a  well-mixed  atmosphere  the  removal  rate,  A  ,  is  determined  by 

R 


(2.7) 


where,  A  =  Removal  rate  (sec  ) 

R 

v  =  Deposition  velocity  (cm/sec) 
d 

2 

A  =  Containment  surface  area  (cm  ) 
s 

3 

V  =  Containment  volume  (cm  ) 


The  deposition  velocity,  v  ,  can  be  analytically  determined 

using  Stoke's  law[15].In  addition,  numerous  experiments  have  been 

conducted  to  determine  v  (and  A_)  for  various  atmospheres,  concentra- 

d  R 

tions,  and  surfaces [6,17] .These  experiments  have  confirmed  that  the 
depletion  of  iodine  follows  the  exponential  relation 

-v 

C ( t)  =  Coe  K  (2.8) 

3 

where  C(t)  =  Concentration  at  time  t  (Ci/cm  ) 

3 

C0  =  Initial  concentration  (Ci/cm  ) 

A  =  Removal  rate  (sec  ^) 

K 

t  =  Time  (sec) 

Depletion  continues  until  the  equilibrium  value  of  1%  of  the  original 
concentration  is  reached. 

The  same  concept  can  be  extended  to  particulate  fission  products. 
These  deplete  primarily  through  agglomeration  into  0.05  to  15  micron 


diameter  particles  and  subsequent  deposition  due  to  gravitational 
settlement.  Particulates  can  also  be  entrained  with  and  condensed 
in  the  steam,  however  this  is  expected  to  have  only  a  slight  effect 
in  the  MITR,  since  the  amount  of  vapor  formed  will  be  small. 

Assuming  that  the  halogens  are  present  in  containment  as  a 
gas  and  all  other  fission  products  as  solids,  the  effectiveness  of 
these  removal  methods  for  the  MITR  can  be  estimated  as  follows. 

Although  the  deposition  velocity  for  iodine  is  dependent  on  concentra¬ 
tion,  temperature  and  surface  type,  experiments  show  that  the  v.  for 

d 

-1  -3 

painted  surfaces  range  from  1.8  x  10  to  3.1  x  10  cm/sec [6,17] . 

From  Appendix  B,  the  surface  area  to  volume  ratio  for  the  MITR 
containment  is  '■'•'0.1.  Using  Eq.  2.7,  this  sets  the  limits  on  the 
removal  rate  from  2.13  to  0.037  hr  respectively.  From  Eq.  2.8, 
this  predicts  two-hour  concentrations  of  1.4%  and  92.9%  of  the 
initial  concentration.  WASH-1400 [15] suggests  the  use  of  X  =  1.38  hr  * 
for  iodine  which  results  in  a  two-hour  level  of  6.3%.  Other  data  in 
the  literature [6,17]  is  reported  in  terms  of  half-lives,  ranging  from 
10  to  30  minutes.  These  result  in  two-hour  levels  of  1%  (reaches 
equilibrium)  and  6.25%.  Accepting  WASH-1400's  value  and  levelizing  over 
the  two-hour  period  yields  an  average  release  fraction  of  33.9%  McCauley[5] 
suggests  the  use  of  25%.  A  value  of  30%  will  be  used  here.  This  is 
expected  to  be  conservative  since  the  natural  circulation  in  the  MITR 
is  expected  to  be  greater  than  that  from  which  the  WASH-1400  values 
were  developed,  leading  to  increased  deposition. 


The  release  fraction  for  solids  can  be  obtained  in  a  similar 

manner  using  the  WASH-1400  particulate  value  of  A  =  0.13  hr  ^[15] 

R  • 

This  predicts  a  two-hour  concentration  of  77.1%  and  a  levelized  value 
of  88.0%.  Since  gravitational  deposition  is  a  slow  process  and  the 
MITR  containment  is  relatively  small,  there  will  be  less  opportunity 
for  particles  to  agglomerate  and  settle  before  reaching  the  walls 
and  possible  leakage  points.  Therefore  a  conservative  value  of  90% 
will  be  used. 

2.2.5  Summary  of  Release  Fractions 

Multiplication  of  the  individual  release  fractions  F  ,  F_,  F  and 

s  '  f  p 

F^  yields  the  overall  fraction  of  the  total  core  fission  product 
inventory  released  to  and  remaining  available  in  the  containment 
atmosphere.  Individual  and  resulting  total  release  fractions,  F^  ,  for 
the  various  categories  of  isotopes  are  summarized  in  Table  2.2. 

It  should  be  noted  that  the  resulting  release  fractions  are  not 
in  keeping  with  current  regulations [ 9 , 10 , 11 ]  .  However,  the  NRC  has 
recognized  [16 , 18 ]that  the  TID  source  term  is  overly  conservative, 
especially  with  regards  to  the  release  of  iodine,  and  several 
studies  are  currently  underway  which  should  result  in  substantiating 
lower  release  fractions  more  in  line  with  this  study  and  the  observed 
releases  from  historical  accidents.  Until  such  results  are  available, 
Stone  and  Webster  Engineering  Corporation  (SWEC)[19]has  suggested 
a  conservative  interim  source  term  based  on  preliminary  results.  Table 
2.2  presents  this  release  alongside  the  present  and  past  MITR  release 


TABLE  2.2  (Cont. 
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fractions.  As  can  be  seen,  the  release  fractions  proposed  in  this 
study  are  more  conservative  than  the  SWEC  release,  except  for  the 
less  important  "solids".  Since  the  true  source  term  is  expected  to 
be  smaller  than  the  SWEC  value,  the  release  obtained  in  this  thesis 
is  expected  to  be  a  definite  overestimate  of  the  actual  release. 


2.2.6  Comparison  with  Historical  Accidents 

Four  accidents  involving  melting  of  U-Al  fuel  plates  in  light 
water  due  to  flow  blockage  have  occurred.  The  Westin^iouse  Test 
Reactor  partial  melting  of  one  element  on  3  April  1960  released  5000 
curies  of  fission  products  from  the  fuel,  260  curies  of  Xe  and  Kr 
immediately  to  the  containment,  and  800  curies  of  activity  in  all  [6,20]. 
The  Engineering  Test  Reactor  had  18  plates  in  6  elements  partially 
melt  on  12  December  1961.  This  released  42  curies  to  the  coolant, 

6  curies  of  gas  and  0.4  of  particulate  to  the  stack [6, 20, 21  ]  .  Releases 
also  occurred  in  the  Materials  Test  Reactor,  13  November  1962,  and  the 
Oak  Ridge  Research  Reactor,  1  July  1963,  which  were  not  as  well 
instrumented  but  which  in  neither  case  resulted  in  a  hazardous 
radiological  release  [6,20]. 

The  other  accidents  of  interest  are  the  SL-1  and  TMI-2  accidents. 

The  Stationary  Low  Power  Reactor  No.  1  underwent  a  reactivity  excursion 

on  3  January  1961,  which  resulted  in  20%  core  melt  with  up  to  10%  of 

the  core  fission  product  inventory  being  ejected  outside  of  the  vessel. 

However,  only  0.01%  of  the  fission  products  escaped  the  reactor  building 

131 

and  only  0.5%  of  the  I  [18,20].  The  Three  Mile  Island  Unit  2 


experienced  a  well-documented  accident  on  28  March  1979  which  resulted 


in  extensive  core  damage.  The  releases  of  fission  products  from  the 
fuel  are  estimated  to  be  50-70%  of  the  noble  gases,  40%  of  the  iodine 
and  2-3%  of  the  solids,  with  releases  to  the  atmosphere  of  10%  for  noble 
gases,  0.00002%  for  iodine  and  a  negligible  amount  of  solids. 

The  amount  of  iodine  retained  in  the  primary  loop  was  500,000  times 
the  amount  released  to  the  atmosphere  [18]. 

From  the  historical  data  it  is  believed  that  the  release  fractions 
in  Table  2.2  are  indeed  conservative. 

2.3  Comparison  With  the  SAR 

Besides  the  inclusion  of  additional  fission  products,  this 

study  differs  from  the  SAR  in  that  the  release  fraction  for  the  noble 

gases  has  been  doubled  and  that  of  iodine  decreased  slightly.  The  net 

14 

effect  is  an  SAR  total  fission  product  release  strength  of  1.70  x  10 

14 

Mev/sec  vs.  a  release  strength  of  2.60  x  10  Mev/sec  in  this  study. 


CHAPTER  3 


ATMOSPHERIC  RELEASE 

3.1  Release  Types  and  General  Assumptions 

Atmospheric  releases  from  reactor  containments  fall  into  two 
categories;  elevated  releases,  typically  from  a  stack,  and  non¬ 
elevated  releases,  usually  due  to  leakage  through  small  cracks  or 
other  penetrations  in  the  containment  building  itself.  During  normal 
operation  of  the  MITR-II,  the  off-gas  treatment  system  results  in 
small  releases  through  the  exhaust  ventilation  and  out  the  stack. 

During  an  accident  the  exhaust  plenum  to  the  stack  will  be  automatically 
closed  and  sealed.  With  the  containment  sealed,  a  rise  in  contain¬ 
ment  pressure  is  possible  depending  on  the  nature  of  the  accident. 

The  design  limit  of  the  containment  is  2  psi  above  atmospheric 
pressure.  If  the  containment  pressure  should  approach  this,  the 
pressure  can  be  reduced  through  use  of  the  pressure  relief  system. 

This  system  consists  of  a  set  of  valves  and  filters  which  will  allow 
controlled  venting  of  containment  atmosphere  through  the  stack. 

If  utilized,  the  filters  can  be  expected  to  retain  about  99%  of  the 
iodines  and  particulate  matter  but  virtually  none  of  the  noble 
gases  (annual  tests  show  a  retention  of  99.8  %)  [1], 

As  the  building  pressure  rises,  leakage  through  small  cracks 
in  the  joints  and  seals  can  occur.  The  leakage  rate  is  proportional 
to  the  amount  of  over-pressure  with  the  maximum  leakage  rate  set 
in  the  technical  specifications  as  1%  of  the  containment  volume  per 
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pound  overpressure  per  day  [1]. 

It  is  not  anticipated  that  the  design  basis  accident  or  any 
other  credible  accident  could  produce  a  rise  in  containment  pressure 
of  2  psi.  For  the  DBA  as  postulated,  and  including  a  slight  decrease 
in  atmospheric  pressure,  a  best  guess  estimate  of  the  pressure  rise 
is  on  the  order  of  0.5  psi.  Therefore  use  of  the  pressure  relief 
system  is  not  expected  for  the  duration  of  the  accident  and  no 
atmospheric  releases  from  the  stack  will  be  considered.  It  is 
assumed  that  some  release  will  occur  due  to  building  outleakage, 
and  it  is  this  release  which  will  be  considered  in  the  following 
sections. 

3.2  Leakage  Rate 

The  estimated  leakage  rate  from  containment  is  dependent  on 
three  conservative  assumptions.  First,  the  containment  pressure 
is  assumed  to  undergo  an  instantaneous  increase  to  the  design  limit 
of  2  psi  over  atmospheric.  As  stated  previously,  this  is  not 
believed  likely  to  occur.  Second,  the  maximum  leakage  rate  of  1% 
of  the  building  volume  per  psi  over  pressure  per  day  will  be  assumed. 
Leakage  tests  of  the  containment  have  shown  actual  leakage  rates 
to  be  lower  [3].  Third,  the  leakage  rate  will  be  assumed  to  remain 
constant  for  the  duration  of  the  accident  when  in  actuality  it  would 
decrease  slightly.  With  the  above  conservative  assumption,  the  leakage 
rate,  X^,  is 
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A  =  . 02V/day  =  2.3  x  10  V/sec 


where  V  =  Volume  of  containment. 


3.3  Leakage  Diffusion  Models 


The  concentration  of  radioactivity  downwind  from  the  contain¬ 
ment  due  to  atmospheric  dispersion  can  be  estimated  from  [22]: 


X/Q  = 


UTO  O  +  cA  U 

y  z  xs 


sec/m" 


(3.1) 


where 


X/Q  =  Relative  concentration  (sec/m  ) 

3 

X  =  Concentration  of  radioactivity  (Ci/m  ) 

Q  =  Release  rate  of  fission  products  (Ci/sec) 
u  =  Wind  speed  (m/sec) 

Oy  =  Lateral  plume  dispersion  coefficient  (m) 

Oz  =  Vertical  plume  dispersion  coefficient  (m) 
c  =  Shape  factor 

2 

A.  =  Containment  cross-sectional  area  (m  ) 
xs 

As  seen  in  Eq.  3.1,  the  initial  concentration  is  reduced  by  two 

contributions,  atmospheric  diffusion  and  the  building  wake  effect. 

The  diffusion  term  (utto  a  )  reflects  the  decrease  in  concentration  due 

y  z 

to  plume  spread  horizontally  and  vertically.  This  is  a  function  of 
the  wind  speed,  distance  downwind  and  atmospheric  stability  conditions. 
Stability  conditions  are  usually  expressed  using  the  Pasquill  stability 
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categories:  A(very  unstable)  through  G  (very  stable).  Value? 
of  a  and  0  as  a  function  of  distance  downwind  and  stability 

y  z 

category  are  available  in  the  literature  [7,  10,  11,  22,  23). 

The  wake  effect  reflects  the  tendency  for  the  released  particles 
to  quickly  expand  to  fill  a  volume  of  air  on  the  order  of  the 
containment  volume.  This  effect  is  a  function  of  windspeed,  building 
cross-sectional  area  and  a  shape  factor  c.  The  shape  factor 
depends  on  the  particular  situation  but  generally  varies  from 
to  2  [22].  As  a  conservative  measure  the  NRC  requires  power  plant 
calculations  to  use  the  minimum  cross-sectional  area  and  a  shape  factor 
of  h. 

While  Eq.  3.1  provides  the  "exact"  value  of  X/Q,  the  NRC 
states  that  the  reduction  of  X/Q  due  to  the  wake  effect  can  be  no 
more  than  a  factor  of  3,  so  that  the  result  obtained  using  Eq.  3.1 
must  be  compared  with  that  from 


X/Q 


1 

3u7T0  o 

y  z 


(3.2) 


and  the  greater  value  used  in  consequence  evaluations  [10,11], 

In  addition,  the  NRC  recently  concluded  that  the  above  equations 
result  in  concentration  estimates  that  are  too  high  during  light 
wind  and  stable  or  neutral  atmospheric  conditions  [23].  For  these 
conditions  lateral  plume  meander  is  important  an'-’  .  an  be  considered 
using 

X/Q  =  - i - 

uMtto  a 

y  z 


(3.3) 
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where  M  =  Meander  correction  factor 

The  correction  factor  M  is  a  function  of  windspeed  and  stability 
category  and  is  shown  in  Fig.  3.1. 

For  the  conditions  above  the  greater  value  of  X/Q  obtained 
using  Equations  3.1  and  3.2  is  compared  to  that  obtained  using 
Eq.  3.3  and  the  lesser  value  selected  [23]. 

3.4  Meteorological  Data 

In  order  to  use  the  concentration  equations  the  site  wind- 
speed  and  stability  category  must  be  known.  The  best  windpspeed 
data  presently  available  is  from  the  National  Weather  Service 
office  at  Boston's  General  Logan  Airport.  The  number  of  recorded 
instances  of  windspeeds  and  their  directions  for  the  year  1981  are 
listed  in  Appendix  A,  Table  A. 2.  The  frequency  of  stability  conditions 
can  be  obtained  from  historical  data  in  the  MITR  Final  Hazards 
Report  [2].  This  information  has  been  summarized  in  Table  A. 3. 

The  appropriate  windspeed  to  use  in  the  calculations  was  deter¬ 
mined  for  each  direction  by  selecting  the  lowest  windspeed  in  that 
direction  which  was  reached  or  exceeded  99.5%  of  the  total  time  (for 
all  directions).  Resulting  directional  windspeeds  are  listed  in 
Table  3.1. 

The  stability  condition  determines  which  graphs  will  be  used  to 

obtain  values  for  c  ,  o  and  M.  The  value  of  X/Q  increases  as 

y  z 

stability  increases.  Therefore  the  stability  condition  which  was  equal 


Stability 

Class 


Figure  3.1 


Correction  Factors  for  Pasqui 11-Gi f ford 
0  Values  by  Atmospheric  Stability  Class  [23] 
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Co  or  less  stable  than  the  most  stable  condition  for  99%  of  the  time 
has  been  selected.  This  corresponds  to  Pasquill  Category  F  (moderately 
stable) . 

3.5  Exclusion  Area  Distances 

In  order  to  determine  0  and  a  the  distance  from  the  contain- 

y  z 

ment  shell  to  the  exclusion  area  boundary  for  each  wind  direction 

must  be  known.  These  values  were  taken  from  an  MIT  Physical  Plant 

site  blueprint  [24]  on  which  the  area  around  the  reactor  was  divided 

into  16  sectors  of  45°  each,  centered  on  each  wind  direction.  The 

shortest  distance  between  the  reactor  containment  shell  and  the 

exclusion  area  boundary  within  each  sector  has  been  designated  the 

sector  distance,  X.  These  values  are  listed  in  Table  3.1. 

In  addition,  the  cross-sectional  area,  A  ,  of  the  containment 

xs 

2 

building  (from  Appendix  B)  is  314  m  . 

3.6  Dispersion  Coefficients 

The  dispersion  coefficients  O  and  a  are  somewhat  difficult 

y  z 

to  determine  as  they  are  given  in  the  literature  as  a  set  of  curves 
2  5 

over  the  range  10  $  x  5  10  meters  and  it  is  impossible  to  accurately 

extrapolate  them  down  to  the  range  of  interest,  8  S  x  $  25  meters. 

One  alternative  is  to  use  the  interpolation  formulas  for  cr  and  a 

y  2 

developed  by  Briggs  which  fit  the  Pasquill  curves.  For  type  F  stability 
they  are  [ 7 ] : 
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Oy  =  0.04  x  (1  +  O.OOOlx)"^  m 

O  =  0.16  X  (1  +  0 . 0003x)_1  m 

where  x  =  Distance  downwind  (m  ) 


(3.4) 

(3.5) 


Since  these  functions  were  fitted  to  the  Pasquill  curves  they 

2  5 

also  are  strictly  valid  only  in  the  range  10  -  x  -  10  meters. 

The  validity  of  the  Briggs  equations  for  shorter  distances  was 

investigated  by  comparison  with  the  older  Sutton  diffusion  parameters 

C  and  C  using  the  conversion  relationships  [22] 
y  z 


a  C?  I'  f 
y  '  7T  x  * 


m 


Cz  1-  ~ 

z  =  7T~  x 


m 


The  difficulty  with  using  Sutton's  equations  is  that  there  is 
no  agreement  among  researchers  as  to  the  correct  values  of  and 
for  a  given  atmospheric  condition.  One  representative  set  of  values 
for  type  F  conditions  are  [25]: 

C  =  0.10 

y 

Cz  =  0.06 
n  =  0.5 


Values  of  a  and  a  calculated  with  the  Briggs  equations  and 

y  2 

the  Sutton  equations  within  the  range  of  interest  are  very  close. 
Therefore  the  Briggs  equations  have  been  used  under  the  assumption 


that  they  yield  the  most  accurate  values  within  the  available  research 
in  this  area  and  that  the  error  involved  is  on  the  same  order  as 


that  of  using  airport  wind  data.  Values  of  a  and  a  evaluated 

y  z 

using  Eqs.  3.4  and  3.5  are  listed  in  Table  3.1. 

3.7  Application  of  the  Diffusion  Models 

Substitution  into  Eq.  3.1  shows  that  the  wake  effect  is  dominant 

over  the  diffusion  term,  so  much  so  as  to  make  the  dispersion  due 

to  diffusion  almost  negligible.  Consequently,  Eq.  3.1  will  always 

predict  a  lower  X/Q  value  than  Eqs.  3.2  and  3.3  which  are  diffusion 

oriented.  Fig.  3.1  gives  the  value  of  M  in  Eq.  3.3  as  4  for' all 

directions.  Therefore,  for  a  given  set  of  values  for  u,  a  and  a 

y  z 

Eq.  3.3  will  yield  a  lower  value  of  X/Q  than  Eq.  3.2  and  should 
be  selected  to  evaluate  X/Q  for  each  sector  as  specified  in  Sec.  3.3. 

While  selection  of  Eq.  3.2  is  in  line  with  the  NRC  guidelines 
it  is  possible  that  this  will  yield  a  X/Q  (and  a  subsequent  dose) 
that  is  too  conservative.  For  example,  for  the  south  sector,  Eq.  3.2 
predicts  a  X/q  0f  1.76  while  Eq.  3.1  with  the  full  wake  effect 
yields  a  X/Q  of  0.00573,  a  factor  of  n,  300  difference.  Due  to  the 
channeling  effect  of  the  adjacent  buildings  it  is  expected  that  the 
turbulent  mixing  of  the  wake  effect  for  the  MITR  will  be  more  pronounced 
than  that  of  a  containment  located  in  open  countryside.  Therefore 
X/Q  will  be  calculated  using  both  the  ’’exact"  Eq .  3.1  and  the  "con¬ 
servative"  Eq.  3.2,  with  the  actual  value  of  X/Q  for  the  MITR  expected 


to  be  closer  to  the  "exact"  value.  Resulting  X/Q  values  for  each 
sector  are  shown  in  Table  3.1. 

The  largest  value  of  X/Q  is  the  controlling  value  for  evaluation  of 
the  leakage  dose  received  at  the  boundary.  As  seen  in  Table  3.1 
this  occurs  at  a  distance  of  8  meters  in  the  direction  south,  which 
is  the  rear  fence  of  the  reactor  site.  This  fence  faces  a  little- 
used  railroad  siding  which  is  seldom  occupied  by  the  public.  The 
boundary  which  is  most  frequently  occupied  by  the  public  is  the 
fence  bordering  the  Albany  Street  sidewalk.  This  sidewalk  has  an 
almost  constant  flow  of  pedestrians.  The  largest  X/Q  value  along  this 
boundary  corresponds  to  the  direction  north  with  a  distance  of  21 
meters.  Although  strictly  speaking,  the  rear  fence  is  the  controlling 
boundary,  the  front  fence  dose  is  of  more  practical  value.  There¬ 
fore,  the  dose  received  at  both  8  and  21  meters  will  be  evaluated 
in  this  and  the  following  chapters. 

3.8  Total  Activity  Release 

The  quantity  of  any  fission  product  initially  available  for 
release  from  containment  is  equal  to  the  saturated  core  fission 
product  inventory,  Q^,  times  the  overall  release  fraction  for  that 
isotope,  F* . 

If  the  fission  product  inventory  is  assumed  to  be  distributed 
uniformly,  then  the  concentration  of  fission  products  is  equal  to 
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scavenging. 


9 


If  the  release  occurs  during  rain  or  snow,  precipitation 


scavenging  would  occur  and  wash  out  some  particles  to  the  ground 


a 


Since  this  would  lower  the  inhalation  (thyroid)  dose,  it  will  be 
conservatively  assumed  that  no  precipitation  occurs  for  the  duration 
of  the  release. 

3.10  Beta  Dose 

The  dose  rate  in  air  from  an  infinite  uniform  cloud  of  beta  radiation 
is  determined  from  [22]: 

D'  _  (1.6  x  10  ^ergs/Mev)  (3 . 7  x  10^dis/Ci-s) EpX 
^  (1293  g/m^)(100  erg/g-rad) 

which  reduces  to: 

3°  =  0.457  E.X  rads/sec  (3.8) 

p 

where, 

D ' 

8  =  Beta  dose  rate  (rads/sec) 

E.  =  Average  beta  energy  per  disintegration  (Mev/dis) 
p 

3 

X  =  Concentration  of  beta-emitting  isotope  (Ci/m  ) 

3 

1293  =  Density  of  air  at  S.T.P.  (g/m  ) 

The  surface  body  dose  rate  is  about  one-half  of  this  [22]  or 

6D  =  0.23EoX  rads/sec  (3.9) 

P 

Since  the  concentration  is  time-dependent  the  total  beta  dose  is: 

7  2  0  0 

3d  =  0 . 23E  f  X(t)dt  rads  (3.10) 

p; 


The  concentration  integral  can  be  related  to  the  previously 
determined  X/Q  values  by  the  relationship 


X(t)  =  (X/Q)0(t)  Ci/m 


which  when  integrated  yields 
7200  7200 


/  X(t)dt  =  (X/Q)  /  Q(t)dt 
0  o 

720  0 

Since,  by  definition  Q  =  J Q(t)dt,  the  above  can  be  expressed 
72  0  0  0 

/  X(t)dt  =  (X/Q)Q 


as 


Ci-sec 


m 


The  total  beta  dose  received  due  to  isotope  i  is  therefore 


gD.  =0.23  E^Q*  (X/Q)  rads 


(3.1 


-i  .. 


The  value  of  E_  is  equal  to  one-third  the  value  of  the  maximum  beta 
p 


energy  for  isotope  i  [22].  Values  of  for  each  isotope  are  listed 

p 


in  Table  A.  4. 

These  conservative  calculations  indicate  an  estimated  total 
beta  exposure  due  to  all  isotopes  of  3.42  rads  at  8  meters  and 


0.517  rads  at  21  meters  using  the  "conservative"  X/Q  values  and 


1.11  x  10  rads  at  both  boundaries  using  the  "exact"  X/Q. 


3.11  Gamma  Doses 


Developed  in  the  same  manner  as  the  beta  equation,  the  dose 


rate  due  to  gamma  radiation  in  air  is  [22]: 
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D'  =  0.507  E  X  rads/sec  (3.12) 

Unlike  the  short-range  beta  radiation,  the  gamma  radiation  received 
is  reduced  by  one-half  due  to  the  presence  of  the  ground  [22], 
resulting  in 

D'  =  0.25  EyX  rads/sec  (3.13) 

The  total  gamma  dose  due  to  isotope  i  is 

yD.  =  0.25  Ey  Q^  (X/Q)  rads  (3.14) 

The  value  of  Ey  in  the  above  equation  is  determined  by  weighting 
the  gamma  energy  spectrum  for  each  isotope  i. 

A  second,  more  accurate,  method  to  determine  the  gamma  dose 
is  to  use  WASH-1400  [7]  computer  generated  dose  conversion  factors, 
Cy,  in  the  equation 

yDi  =  CyQ*(X/Q)  rads  (3.15) 

where  Cy  =  Photon  dose-conversion  factor  for  immersion  in  air 
contaminated  with  isotope  i  (rads  per  iL1.~|.e£  ) 

m 

i 

Values  of  Cy  are  not  listed  for  all  isotopes  of  interest. 

Where  unavailable,  values  of  E^were  determined  and  calculated 

y  i 

using  Eq.  3.14.  Values  of  Ey  and  Cy  are  listed  in  Table  A. 4.  These 

conservative  calculations  indicate  an  estimated  total  gamma  dose 

due  to  all  isotopes  of  4.70  rads  at  8  meters  and  0.710  rads  at  21 

-2 

meters  using  the  "conservative"  X/Q  values  and  1.55  x  10  rads  at 
both  distances  using  the  "exact”  X/Q. 
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3.12  Thyroid  Dose 

In  general,  the  radiation  dose  received  by  the  thyroid  is  de¬ 
pendent  on  the  amount  of  isotope  inhaled,  the  fraction  which  is 
transported  to  the  thyroid,  the  energy  absorbed  by  the  thyroid 
per  disintegration,  and  the  rate  at  which  the  isotope  is  removed 
by  decay  and  biological  elimination. 

WASH-1400  [7]  also  provides  inhalation  dose  coversion  factors 
for  various  parts  of  the  body,  including  the  thyroid,  which  can 
be  used  in: 

T°i=  BrCTQT  (X/Q)  rads  (3.16) 

where  jDi=  Dose  to  thyroid  from  isotope  i  (rads) 

3 

B^_=  Breathing  rate  (m  /sec) 

C^=  Thyroid  inhalation  conversion  factor  for  isotope  i 
(rads  per  Ci  inhaled) 

Isotopes  of  interest  which  are  not  inlcuded  in  WASH-1400  were 
checked  against  ICRP  Report  // 2  [26]  and  found  to  have  no  contribution 
to  the  thyroid  dose. 

The  standard  breathing  rate  for  calculations  of  internal  dose 
is  3.47  x  10  4  m^/sec  [23]. 

Values  of  for  each  isotope  are  listed  in  Table  A. 4.  The  result¬ 
ing  "conservative"  estimate  of  the  total  thyroid  dose  due  to  all 
isotopes  is  36.3  rads  at  8  meters  and  5.49  rads  at  21  meters  and  the 
"exact"  estimate  is  1.18  x  10  3  rads  at  either  boundary. 


3.13  Summary  and  Comparison  with  the  SAR 


As  shown  in  Table  3.2,  the  conservative  NRC  atmospheric  release 
calculations  performed  in  this  chapter  indicate  that  the  estimated 
maximum  total  potential  for  personnel  exposure  due  to  containment  leakage 
over  a  two  hour  period  following  the  accident  is  8.12  rads  whole-body 
and  36.3  rads  to  the  thyroid  at  the  back  fence  (  a  distance  of  8 
meters)  and  1.23  rads  whole-body  and  5.49  rads  thyroid  at  the  Albany 
St.  fence  (a  distance  of  21  meters).  The  more  realistic  estimate 

-2 

based  on  the  full  wake  effect  predicts  a  maximum  dose  of  2.66  x  10 
rads  whole-body  and  1.18  x  10  ^  rads  to  the  thyroid  at  either 
boundary.  A  comparison  of  the  "conservative"  and  "exact"  whole-body 
and  thyroid  doses  over  the  entire  range  of  exclusion  area  distances 
is  shown  in  Figs.  3.2  and  3.3.  Note  that  the  wake  effect  is  so  large 
as  to  make  the  difference  in  "exact"  doses  due  to  distance  negligible. 

The  SAR  gives  the  two-hour  leakage  dose  to  the  thyroid  (whole-body 
was  not  evaluated)  as  2.8  x  10  ^  rads  [3].  This  was  determined  using 

X/q  =  — 1 —  sec/m^ 

cA  u 
xs 

in  which  full  credit  is  taken  for  the  wake  effect  and  diffusion  is 
neglected.  In  addition,  the  SAR  assumed  simultaneous  full  operation 
of  the  pressure  relief  system  and  a  windspeed  of  1  m/s. 


TABLE  3.2 

Leakage  Dose  Summary 


Dose 


Component 

of 

Dose 

"Exact" 
Eq.  3.1 

(Rads) 

"Conservative" 
Eq.  3.2 

8m  21m 

Whole-body : 

Beta 

0.0111 

3.42 

0.517 

Gamma 

0.0155 

4.70 

0.710 

Total 

0.0266 

8.12 

1.23 

Thyroid: 

0.118 

36.3 

5.49 

Figure  3.2  Leakage  Whole-body  Doses 


Figure  3.3  Leakage  Thyroid  Doses 
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CHAPTER  4 
DIRECT  GAMMA  DOSE 


4.1  General 


Those  isotopes  which  do  not  leak  from  the  containment  will 
constitute  a  source  of  gamma  radiation.  The  concrete  shadow  shield 
is  designed  to  attenuate  this  radiation  and  reduce  the  hazard  to 
individuals  near  the  containment.  Due  to  the  short  exclusion  area 
distances  involved,  the  gamma  dose  from  radiation  penetrating  the 
shadow  shield  will  be  non-neglible.  In  addition,  the  fission  products 
located  in  that  portion  of  the  containment  volume  above  the  shadow 
shield  are  only  shielded  by  the  steel  shell.  In  the  following 
sections  the  dose  resulting  from  these  two  penetrations  will  be 
determined . 


4.2  Gamma  Source  Term 

The  gamma  source  within  the  containment  will  consist  of  those 
fission  products  which  deposit  or  plateout  and  those  airborne  fission 
products  which  do  not  leak  out.  As  stated  in  Sec.  3.8  the  initial 
quantity  of  fission  product  i  airborne  in  the  containment  is  equal 
to  Q^.  This  will  be  reduced  over  time  due  to  leakage  and  decay. 
The  quantity  which  deposits  inside  the  containment  is  equal  to 

F1F^F1(l-Fi)Q1 
s  f  p  c  xs 


which  can  also  be  expressed  as 


(1-FC>  Fi0i 

Vs 

c 


The  quantity  of  deposited  fission  products  will  decrease  over  time 
only  by  radioactive  decay. 

The  time-dependent  containment  inventory  of  fission  product  i 


is  therefore 


q! (t)  =  F^qV(Xl  +  Ai)t  +  (4-  -  1)  fV>  xic  Ci 


.  .  -(A.  +A.)t  -  -A.t 

Qc(t)  =  FRQse  1 


Ci  (4.1) 


The  total  number  of  decay  emissions  from  isotope  i  over  the  two 
hour  period  is  given  by 

72  0  0 

Q*  =  (3.7  x  10^  dis/sec-Ci)  [  Q*(t)dt 


which  after  integration,  is 


-(A  +  A . ) 7200  -A 

Q*  =  3.7  x  1010  f£QJ  ^  ~T~~TTT -  +  (Ii  "  1)J-=-V 

L  l  F  l 

L  c 


A. 7200 

l 


dis  (4.2 


Resulting  values  of  Qc^  for  each  isotope  are  listed  in  Appendix  A,  Table 


The  energy  and  abundance  of  each  isotope's  gamma  decay  spectrum 
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are  also  shown  in. Table  A. 5.  For  convenience,  photons  have  been 

grouped  into  discrete  energies  following  a  logrithmic  scale,  with 

individual  photons  being  allocated  to  .the  closest  energy. 

The  total  number  of  emissions  for  each  energy  is  then  equal  to 

the  product  of  the  number  of  emissions,  Q1  ,  for  each  isotope  and 

CT 

the  photon  abundance  for  that  isotope  at  that  energy,  summed  over  all 
isotopes.  The  resulting  total  number  of  gammas  for  each  energy  is 
divided  by  the  containment  volume  and  duration  of  release  to  obtain 
the  time-averaged  total  containment  volumetic  source  strength, 

SVT>  Values  of  SvT  for  each  energy  E  are  listed  in  Table  4.1. 


4.3  Direct  Dose  Modeling 


From  any  given  point  on  the  ground  the  containment  can  be 
divided  into  two  volumes,  that  which  lies  above  the  plane  passing 
through  the  top  edge  of  the  shadow  shield  and  the  detector  and  that 
which  lies  below.  These  volumes  will  be  designated  and  respect¬ 
ively,  as  shown  in  Fig.  4.1. 

That  portion  of  the  containment  volume  which  is  contained  in 
can  be  estimated  by  approximating  as  one  half  of  an  ellipsoid  with 
dimensions  a,  d,  and  e  as  shown.  The  volume  may  then  be  obtained 


=  1/2  x  -j  tt  ade 


=  -jfade 


(4.3) 
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TABLE  4.1 

Containment  Volume  Source  Strength 

Gamma  Energy  Total  Containment  Volume 

Source  Strength 
b  S 

V 

(Mev)  T 

(Photons /cm3- sec) 


0.03 

0.04 

0.05 

0.06 

0.08 

0.10 

0.15 

0.20 

0.30 

0.40 

0.50 

0.60 

0.80 

1.0 

1.5 

2.0 

3.0 

4.0 


3.35  x  102 
6.11  x  10~ 
2.16  x  102 
1.96  x  101 
1.25  x  104 
9.90  x  10° 
6.72  x  103 
7.14  x  103 
7.22  x  103 
7.37  x  103 
5.29  x  103 
4.46  x  103 
1.08  x  10  4 

2.69  x  103 

3.70  x  103 
1.04  x  104 
1.18  x  103 
1.32  x  101 
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The  parameters  of  are: 

for  the  back  fence  (8m):  a  =  1.22  m 

d  =  2 .44  m 
e  =  7.62  m 

and  for  the  front  fence  (21  m) :  a  =  1.83  m 

d  =  6.10  m 
e  =  10.7  m 

Substitution  in  Eq.  4.3  yields  the  result: 

V  (8)  =  47.5  m3 

V  (21)  =  250  m3 

3  3 

Since  V  =  4.73  x  10  m  ,  it  is  seen  that 

V  (8)  ^.01  V  (4.4) 

V  (21)  n,.  05V  (4.5) 

Consequently,  since  =  V,  V^can  be  expressed  as 

V  (8)  =  ^.99V  (4.6) 

V  (21)  =  %.95V  (4.7) 

Volume  V  can  be  considered  to  be  a  volumetric  gamma  source 
1 

shielded  by  a  single  steel  slab  and  V^  a  volumetric  source  shielded  by 
a  two  layer  concrete/steel  slab. 

The  easiest  method  available  to  determine  the  dose  from  volume¬ 
tric  sources  is  the  point  kernel  technique.  Since  each  gamma  decay 
and  attenuation  interaction  is  independent  of  any  other,  a  volume 
source  can  be  considered  to  be  a  number  of  point  isotropic  sources. 
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A  ray  is  drawn  from  a  differential  source  point  to  the  detector  and  the 
detector  response  due  to  the  ray  path  determined.  This  is  the  kernel 
response.  The  effect  due  to  all  points  is  obtained  by  integrating 
over  the  volume  of  the  source.  Point  kernel  integrations  have  been 
performed  for  a  number  of  standard  shapes  and  are  available  in  any 
of  the  standard  shielding  manuals  [27,  28,  29]. 

In  this  study  the  upper  volume  will  be  treated  as  a  sphere  of 
volume  located  at  a  point  just  above  the  shadow  shield  and  the 
lower  volume  as  a  right  circular  cylinder  of  volume  with  a 
height  to  radius  ratio  of  one  (see  Fig.  4.2).  Other  source  geometries 
such  as  a  line  source  and  truncated  cone  were  also  considered  but 
found  to  be  less  conservative. 


4.4  Penetration  Calculations 


4.4.1  Steel  Shell  Penetration 


A  spherical  volume  source  of  constant  source  strength  S^ 

3 

(photons/cm  -sec)  can  be  approximated  by  a  disk  of  the  same  radius 
having  a  surface  source  strength 


4  2 

S  =  —  RS  photons/cm  -sec 

A  3  v 


(4.8) 


located  at  a  self-absorption  distance  z  [29]  (see  Fig.  4.3a). 
Assuming  that  the  containment  atmosphere  is  primarily  air,  self¬ 
absorption  will  be  small  and  it  is  conservative  to  assume  z  =0. 

The  volume  then  reduces  to  a  disk  of  radius  located  a  distance  x^ 


from  the  exclusion  area  boundary  point  P  as  shown  in  Fig.  4.3b. 
The  flux  received  at  P  due  to  photons  of  energy  E  from  this 


disk  source  shielded  by  a  steel  slab  is  [29]: 
BS, 


d>  = 

i 


3Al 


W  '  VV600:) 


photons/cm  -sec 


(4.9) 


where : 


<j)^  =  Photon  flux  (photons/cm  -s) 

B  =  Build  up  factor 

SA^  =  surface  source  ^trength  for  volume  and  energy  E 
(photons/cm  -sec) 

01  =  ^ST^ST ^umber  of  mean  free  paths  in  the  steel  shield) 
lJg^,  =  Linear  attenuation  coefficient  for  steel  (cm  ^) 

=  Steel  thickness  (cm) 


CO 


dt 


for  the  derivation  of  Eq.  4.9  see  Appendix  C.  Substituting  for 
from  Eq.  4.8  the  flux  becomes 


=  3  BRlSvl 


E(b^)  -  E^  (b^sec  9) 


photons/cm  -sec 


(4.10) 


Buildup  and  attenuation  in  the  air  will  be  neglected.  Both 
effects  are  small  and  tend  to  cancel  each  other.  Values  of  y  and 
subsequent  values  of  b^  are  shown  in  Tables  A. 6  and  A. 7. 

The  dose  at  P  is  determined  using  the  conversion  factor  in 

Dose  =  rads  (4.11) 
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where 


—  9 

(E  Mev/photon)  (1.6  x  10  ergs/Mev) (Ua  cm  /g) (7200s) 
100  ergs/g-rad 


which  reduces  to 


-4  2 

C.  =  1.15  x  10  Ey  rads/photons/cm  -s 
D  a 


(4.12) 


where  y^  =  True  energy  absorption  coefficient  in  air  (cm  /g) . 

(See  Table  A. 6) 

(Note:  for  convenience,  the  energy  dependency  notation  (E)  has  been 
dropped  in  references  to  S  ,  C  ,  dose,  y  ,  etc.  The  reader  should  keep 

V  1/  cl 

in  mind  that  all  dose  equations  are  for  a  given  photon  energy,  not 
the  total  dose) . 

Substituting  C  and  $  into  Eq.  4.11  the  dose  becomes 

y  Y 

Dose  =  7.67  x  10~5  Ey  BR  Sv  E  (b  )  -  E  (b  sec0  )  rads  (4.13) 

el  X  X  XX  XX  X 

The  buildup  factor  B  (a  function  of  energy  and  mean  free  path 
length)  is  included  to  reflect  the  tendency  during  broad  beam  atten¬ 
uation  for  gammas  not  originally  on  a  direct  path  to  the  detector 
to  scatter  towards  the  detector  while  traversing  the  medium.  It  is 
defined  as  the  ratio  of  the  total  detector  response  to  the  uncollided 
(attenuated)  response.  Tabulated  values  of  exposure  (dose)  and  energy 
deposition  build  up  factors  are  available  for  point  isotropic  and 
broad  beam  parallel  sources  in  an  infinite  medium.  Point  isotropic 
dose  build  up  factors  for  steel  are  shown  in  Table  A. 8.  The  use  of 


infinite  medium  factors  is  conservative  since  build  up  factors  for 
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finite  shield  geometry  are  smaller  [27].  Strictly  speaking,  the  build 

up  factor  values  for  point  sources  cannot,  in  the  case  of  extended 

sources,  be  used  as  simple  multiplicative  factors;  they  can  only  be 

combined  with  the  point  kernel  in  the  integrand.  This  is  because 

the  penetration  paths  of  the  primary  radiation  in  the  shield  are 

always  longer  for  extended  sources  than  the  penetration  paths  of  the 

normal  incident  radiation  from  point  sources,  with  a  consequent 

increase  in  the  build  up  factor  [27]. 

One  way  to  combine  the  build  up  factor  in  the  integrand  is  to 

express  it  as  an  analytic  expression  which  is  compatible  with  the 

integration.  One  of  the  most  common  forms  is  the  Taylor  formula  [27]: 
-o^yT  -a2yT 

B  =  A  e  +  (1-A)  e  (4.14) 

in  which  A,  at ,  and  a2  are  energy-dependant  coefficients  which  are 

fitted  to  the  point  build  up  factor  values.  Substituting  the  ex¬ 
pression  for  B  into  the  equation  for  E^(b)  and  integrating  the  result  is 
E1(b1)  =  AE1(b1’)  +  (1-A)  E1(b1") 

where  b^ '  =  (1  +  a^)b^ 

b1"  =  (1  +  o2)b1 

Therefore,  Eq.  4.13  can  be  rewritten  as 

Dose  =  7.67  x  10~5Ey  R. Sv,  AE  (b  ' )  +  (l-A)E, (b  ") 

all  11  11 

-  AE^(b^'sec8^)-(l-A)E^(b^"sec8^)  rads  (4.15) 

The  function  E  (b)  cannot  be  evaluated  analytically.  Numerically 


E 


determined  graphs  are  available  in  the  literature  but  interpolation 
is  difficult  and  inaccurate.  Fortunately  it  is  possible  to  determine 
the  difference  between  two  E^(b)  functions  using  the  relation  [27] 

El(b)  -  Ei  [bd  +  6)]  =  6e"b  (4.16) 

for  6«1  and  b  >  0. 

For  this  situation,  (1  +  6)  =  secB 

or  6  =  sec0  -  1 

Substitution  into  Eq.  4.15  yields  the  final  result, 

-5  r  -v 

Dose  =  7.67  x  10  Ey  R, SV1  A(sec0,-l)e  +  (1-A)(sec0,  -l)e 

all  1  1 

rads  (4.17) 


where  b^'  =  (1  +  a^)b^ 

b1"  -  (1  +  a2)b1 

Assuming  that  the  fission  products  are  uniformly  distributed  in  the 
containment,  the  volume  relations  in  Eqs.  4.4  and  4.5  lead  to  the 
source  strength  relations 

SV1(8)  =  .01  SVT 

SVl(21)  -  .05  SVT 

The  scattering  geometry  parameters  are: 

0^  =  0.179  radians 
R^  =  2.25  x  lO^cm 


8  m: 
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21  m:  0^  =  0.169  radians 

=  3.90  x  lO^cm 

Values  of  the  coefficients  A,  a^,  and  a ^  are  listed  in  Table 
A. 9  with  the  resulting  doses  listed  in  Table  4.2.  Since  the  Taylor 
coefficients  are  not  available  for  E  <  0.5  Mev,  the  lower  energy 
doses  have  been  computed  using  the  appropriate  point  build  up 
factor  data  from  Table  A. 7  in 

-5  ~k]_ 

Dose  =  7.67  x  10  Ey^R^Sy^B (sec0^  -  l)e  rads  (4.18) 

Doses  for  E  <  .10  have  not  been  determined  as  build  up  factor 

data  for  steel  in  this  energy  range  is  not  available  and  the  increasing 

attenuation  at  lower  energies  makes  the  dose  at  these  energies  negligible. 

As  seen  in  Table  4.2  the  total  dose  at  21  meters  is  larger 

than  that  at  8  meters.  This  is  due  to  the  fact  that  the  portion 

of  the  containment  volume  seen  above  the  shadow  shield  increases 

as  point  P  moves  away  from  the  containment.  An  estimate  of  the 

behavior  of  the  total  dose  byeond  21  meters  was  made  by  tracking 

the  behavior  of  the  largest  dose  component,  that  of  E  =  2.0  Mev. 

The  results,  shown  in  Fig.  4.5,  indicate  that  the  total  steel  pene- 

-2 

tration  dose  goes  through  a  maximum  of  %6.5  x  10  rads  at  a  distance 
of  n,  37  meters.  Beyond  this  the  distance  effect  overcomes  the  increase 
in  volume  and  the  dose  decreases. 

4.4.2 _ Shadow  Shield  Penetration 


As  stated  previously,  that  portion  of  the  containment  volume 


Gamma 


E 

(Mev) 

0.10 

0.15 

0.20 

0.30 

0.40 

0.50 

0.60 

0.80 

1.0 

1.5 

2.0 

3.0 

4.0 


Total 
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TABLE  4.2 

Dome  Penetration  Doses 


Dose 
( Rads ) 

8m  21m 


1.31 

X 

10-8 

1.05 

X 

io-7 

3.80 

X 

10“  5 

2.95 

X 

10“  4 

7.78 

X 

10“5 

6.05 

X 

10“4 

1.46 

X 

10“4 

1.13 

X 

io-3 

2.16 

X 

io“4 

1.68 

X 

io“3 

1.95 

X 

10~4 

1.52 

X 

10“3 

1.99 

X 

io“4 

1.55 

X 

10“  3 

6.36 

X 

10-4 

4.94 

X 

10“3' 

1.92 

X 

io“4 

1.49 

X 

10“  3 

3.62 

X 

10“4 

2.82 

X 

io“3 

1.24 

X 

io"3 

9.63 

X 

10“3 

1.85 

X 

10"  4 

1.44 

X 

10“3 

2.52 

X 

io-6 

1.96 

X 

10“5 

3.49 

X 

io-3 

2.71 

X 

10“  2 

73 


behind  the  shadow  shield  can  be  treated  as  a  right  circular  cylinder 
of  radius  and  height  h^  shielded  by  a  slab  shield  of  thickness  b^ 
as  shown  in  Fig.  4.4. 

For  this  situation  the  flux  at  point  P  is  given  by  [27]: 

BR  Sv  2 

\  -  G  (k,p>ygR2,  b2)  photons/cm  -s  (4.19) 

^2 

where :  k=  — 

R2 


p  =— f—  (must  be  >  1.25) 


y„  =  Linear  attenuation  coefficent  in  the  source  medium 

(cm-1) 

b_  =  y  T  +  y  T  _  =  Total  shadow  shield  thickness  in 
zccblai  -  , 

mean  free  paths 

G  =  Attenuation  function 

The  derivation  of  Eq.  4.19  and  details  of  the  attenuation 
function  G  are  contained  in  Appendix  D.  The  above  equation  can  also 
be  used  for  a  source  with  no  self-absorption  by  setting  ygR2  = 

As  before,  the  dose  at  P  is  found  by  use  of  the  conversion  factor 
Cp,  with  the  result 


Dose  = 


1.15  x  10 


-4 


2  Ti 


Ey a BR? Sy9 


G(k,p,ysR2,b2)  rads 


(4.20) 


The  exact  build  up  factor  for  a  multi  layer  shield  can  only 
be  determined  by  the  solution  of  the  Boltzmann  transport  equation, 
with  appropriate  boundary  conditions,  by  a  numerical  method,  such  as 
Monte  Carlo  technique.  For  most  practical  purposes  simpler  approximate 


Figure  4.4 

Geometry  for  Cylindrical  Volume  Source  with  Slab  Shield 
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methods  are  adequate.  One  method  is  to  use  the  build  up  factor 
data  for  the  last  material  penetrated  based  on  the  thickness  (in 
mean  free  paths)  of  the  entire  shield.  This  method  can  only  be  used 
if  the  last  layer  is  at  least  three  mean  free  paths  thick  [30]. 
Alternatively,  if  almost  all  of  the  attenuation  (measured  in  mean 
free  paths)  occurs  in  a  particular  layer,  then  the  build  up  factor 
is  based  on  that  material  and  the  total  thickness  (in  MFP's)  of  the 
shield  [27] . 

Referring  to  Table  A. 7  one  can  see  that  build  up  factor  data 
for  concrete  should  be  used  for  all  energies  equal  to  or  greater 
than  0.10  Mev.  and  that  for  steel  below  0.10  Mev.  Note  that  the 
concrete  portion  of  the  shadow  shield  is  treated  as  a  layer  of 
pure  concrete  when  in  actuality  it  contains  numerous  steel  reinforc¬ 
ing  bars.  This  is  conservative  as  the  attenuation  in  concrete  is 
less  than  that  in  steel  while  the  build  up  factor  in  concrete  is 
greater  than  in  steel. 

As  before,  greater  accuracy  can  be  achieved  by  incorporating  an 
analytic  build  up  factor  expression  into  the  integrand,  in  this  case 
the  G  function.  Using  the  Taylor  formula  the  dose  now  becomes 
1.15  x  10_4Eu  R?Sv 

Dose  =  - ^ - 3  -  AG(k,p,psR2,b2') 

+  (l-A)G(k,p,UsR2,b2")  rads  (4.21) 

where 

b2'  -  (1  +  ax)b2 


b2"  =  (1  +  a2)b2 


pi 

f- 


H  ■  5  •  .  |<  if  . 
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Values  of  the  Taylor  coefficients  for  concrete  are  listed  in  Table 
A. 9. 

is  related  to  the  total  source  strength  Sv^,  through  Eqs. 

4 . 6  and  4.7. 

The  parameters  of  the  G  function  are  determined  as  follows: 
From  Eqs.  4.6  and  4.7  the  respective  volumes  are 

V2(8)  =  4.68  x  103  m3 

V2(21)  =  4.49  x  103  m3 

For  convenience,  k  is  set  equal  to  one.  This  eliminates  one  set 
of  interpolations  in  the  G  function  tables  and  is  not  too  far  from 
the  actual  containment  h/R  ratio.  Given  that  k  =  1,  and  therefore 
R2  =  h^,  the  radii  can  be  solved  for  using 

V2  =  frR2h2  t0  yield 

R2(8)  =  11.4  m 
R2(21)  =  11.3  m 

Since  s  is  the  total  distance  from  the  center  of  to  P  and  the 
thickness  of  the  shadow  shield  is  0.61  m  (2  ft)  the  variable  p  can 
be  determined  to  be 


P(8)  = 


11.4  +  0.61  +  8 
11.4 


or 


p (8)  =  1.75 
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Since  self-absorption  is  neglected  u  R„  =  0.  Values  of  b  *  and 

s  2  2 

b 2"  are  listed  in  Table  A. 10  along  with  the  corresponding  G  function 
values . 

The  resulting  doses  are  included  in  Table  4.3.  Again,  doses 
for  E  <.10  Mev  have  not  been  determined  for  the  same  reasons  as 
cited  in  Sec.  4.4.1.  As  opposed  to  the  steel  penetration  dose,  the 
shadow  shield  dose  is  a  continuously  decreasing  function  of  distance, 
since  the  volume  decreases  with  distance.  An  estimate  of  the  shadow 
shield  penetration  dose  behavior,  based  on  E  =  2.0  Mev,  is  shown 
in  Fig.  4.5 

4.5  Summary  and  Comparison  with  the  SAR 

The  point-kernel  integration  techniques  used  in  this  section 
to  evaluate  the  direct  gamma  dose  reaching  the  exlusion  area  boundary 

indicate  that  the  estimated  two-hour  dose  due  to  both  steel  and 

-2 

shadow  shield  penetration  is  4.68  x  10  rads  at  the  back  fence 

-2 

(8  meters  from  the  containment)  and  4.76  x  10  rads  at  the  Albany  St. 

fence  (21  meters  from  containment).  The  increasing  steel  penetration 

dose  is  offset  by  the  continuously  decreasing  shadow  shield  dose  to 

produce  the  total  dose  curve  shown  in  Fig.  4.5  with  a  maximum  of 
-2 

^7.6  x  10  rads  at  a  distance  ofn/35  meters. 

-2 

The  SAR  gives  the  total  dose  at  8  meter.,  as  'v  3.1  x  10  rads 

-2 

and  that  at  21  meters  as  'v  2.5  x  10  rads.  The  differences  are  due 
to  the  slightly  different  penetration  methodology  used  in  the  SAR 


as  well  as  the  different  source  terms. 
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TABLE  4.3 

Shadow  Shield  Penetration  Doses 


Gamma  Energy 
E 

(Mev) 


Dose 

(Rads) 


0.10 

1.06 

X 

10 

5.96 

X 

10 

0.15 

6.66 

X 

10"  9 

4.91 

X 

10 

0.20 

1.49 

X 

io“7 

1.09 

X 

10 

0.30 

2.78 

X 

io“6 

1.95 

X 

10 

0.40 

1.56 

X 

io“5 

1.08 

X 

10 

0.50 

3.76 

X 

io“5 

2.63 

X 

10 

0.60 

6.10 

X 

10“  5 

3.55 

X 

10 

0.80 

6.15 

X 

10“4 

4.13 

X 

10 

1.0 

5.23 

X 

10~4 

•  1.98 

X 

10 

1.5 

3.59 

X 

10“  3 

1.70 

X 

10 

2.0 

2.60 

X 

io“2 

1.27 

X 

10 

3.0 

1.23 

X 

io“2 

5.33 

X 

10 

4.0 

2.04 

X 

10“4 

1.11 

X 

10' 

Total 

4.33 

X 

10“2 

2.05 

X 

10' 

ontainment  (M) 


CHAPTER  5 


SCATTERED  GAMMA  DOSE 


5.1  General 

Although  the  shadow  shield  is  quite  effective  in  stopping  direct 
radiation,  the  open  top  construction  leads  to  the  possibility  of 
gamma  radiation  escaping  upwards  and  being  subsequently  scattered 
back  towards,  the  ground  through  interaction  with  the  air  or  the  steel 
dome.  This  scattered  radiation  is  commonly  called  skyshine.  In  the 
following  sections  a  skyshine  analytical  model  will  be  developed 
and  applied  to  scattering  from  the  air  and  the  steel  roof. 


5.2  Scattering  Model 


Consider  a  point  source  of  strength  S  photon/sec  located  a 
distance  x  from  a  detector  at  point  P  with  a  small  scattering  volume 
dV  located  at  a  distance  r^  from  the  source  and  a  distance  r^  from 
the  detector  as  shown  in  Fig.  5.1. 

The  scattering  angle  is  6  =  +  (p  and  the  differential  dose 

recieved  at  P  is 


dD  = 


%SN 


da 


/  2  2 
4l,rl  r2 


dfi 


(0,E)dV 


(5.1) 


where  N  =  Density  of  electrons  in  air  = 


20  ■» 

3.6  x  10  ^lec/cm  at  STP 


X 


Figure  5.1 
Scattering  Geometry 
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^  (9,E)  =  Klein-Nishina  differential  collision  cross- 

section  for  angle  0  and  photon  energy  E. 

2 

(cm  /steradian). 


Eq.  5.1  neglects  build  up  and  attenuation  in  the  air.  There  is 
cylindrical  symmetry  about  x  so  that  [31] 

dV  =  (2ir^simjj)  (r^di^)  (dr^) 

From  the  law’ of  sines  it  is  possible  to  state 


sin<J>  simjj  sin  [it  -  (l[rt-<j>)] 


sin(if>  +  <J>) 


(5.2) 


Differentiating  while  holding  ip  constant  yields 


dr 

d<{> 


i  =  x  sin(ijj  +  cj>) cos -  x  cos(^  +  (}>)  sint]) 
sin^(ijj  +  (j>) 


which  simplifies  to 


dr 


1^  _  x  sinif; 


^  sin^(ij;  +  0) 


Substituting  from  Eq.  5.2,  this  becomes 

2 


dr. 


x  sin  ip 


dip 


dV  can  now  be  written  as 
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dV  =  (2irrlSinW  (r^)  1  d* 


,  2  2 
27Tri  r2 

dV  =  - - — - -  dtJ>d<J> 


(5.3) 


with  the  result  that 


C  SN 

dD  =  2T-  (6.E)dW 


(5. A) 


Dividing  by  a  factor  of  2  to  include  only  the  area  above  the  S-P 
plane  and  integrating  the  total  dose  received  at  P  is 

r  ex,  77 

C  SN  r  r  a 

DoSe  =  ~4x~  j  J  df  rads  l 

oo 

The  Klein-Nishina  differential  collision  (flux)  cross-section 
is  given  by  [32] 

r^  2 

da  _  e  /E*  W  E  E*  .  2a\  2. 

d?2  '  Tlfj  \  "e7"  +  E  "  Sin  Q)cm  /steradian  ( 


(5.5) 


(5.6) 


—  13 

where  r^  =  Classical  radius  of  the  electron  =  2.818  x  10  cm 


E  =  Incident  photon  energy  (Mev) 
E'  =  Scattered  photon  energy  (Mev) 


E  and  E'  are  related  by  [32] 


1  +  — ^Tj-(1  -  cos8) 


(5.7) 
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The  greatest  energy  loss  occurs  when  6  =  180°  and  the  energy  loss  is 
zero  at  9  =  0°.  The  Klein-Nishina  formula  reflects  the  fact  that 
gamma  (Compton)  scattering  becomes  peaked  in  the  forward  direction 
with  increasing  photon  energy. 

5.3  Air  Scattering 

Since  the  higher  energy  photons  will  tend  to  scatter  forward 
with  little  energy  loss  the  containment  volume  will  be  divided  into 
two  source  regions.  The  domed  portion  above  the  shadow  shield,  from 
which  photons  need  only  scatter  through  small  angles  to  reach  the 
ground,  will  be  designated  V  .  The  lower  portion  within  the  shadow 
shield,  which  requires  large  scattering  angles  to  reach  the  ground, 
will  be  designated  V^. 

5.3.1  Upper  Containment 

Assuming  the  fission  products  in  the  dome  are  concentrated  as 
a  point  source  at  the  center  the  scattering  geometry  is  as  shown 
in  Fig.  5.2. 

Eq.  5.5  can  be  applied  directly  if  the  limits  of  integration  are 
changed  so  that  [33] 

TT-<P0  TT-iP 

CnS  N  f  f  Hn 

Dose  =  — ^ - J  dip  J  dip  ^  (9  =\p  +  <£)  rads  (5.8) 
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where  =  Source  strength  in  the  dome  (photons/sec) 

=  Initial  value  of  \p  (radians) 

<pg  =  Initial  value  of  <p  (radians) 

Only  the  area  above  the  S-P  plane  is  considered  as  the  shadow 
shield  will  prevent  interactions  below  this  plane.  In  addition,  the 
above  limits  of  integration  are  conservative  as  the  presence  of 
the  ground  will  prevent  (J>  from  actually  reaching  the  value  TT  -ip. 

No  adjustment  is  made  here  to  reflect  this  since  in  the  next  section 
(lower  containment  scattering)  (p  will  be  able  to  reach  it-  \p. 

Eq.  5.8  can  only  be  evaluated  numerically.  A  simple  method 


is  to  approximate  the  first  integral  by  a  summation  so  that 


:_s  n  v- ■' 

D  u  \ 

4x  zL 


-(6  =  ip  +  tp)  dtp  rads 


(5.9) 


This  allows  the  value  of  i(/  in  0  =  \p  +  <j>  to  be  fixed  and  the 

second  integral  can  be  evaluated  with  a  standard  numerical  technique 

such  as  Simpson's  rule.  The  number  of  intervals  for  \p  and  <p  can  be 

adjusted  to  whatever  accuracy  is  desired. 

The  total  containment  source  strength  S  (photons/sec)  is  calculated 

c 

by  multiplying  the  volumetric  source  strengths  in  Table  4.1  by  the 

volume  of  the  containment.  The  results  are  given  in  Table  5.1. 

Su  can  be  related  to  by  the  fact  that  (from  Appendix  B) 

Vu  ^  0.3V.  If  the  fission  products  are  again  assumed  to  be  distributed 

uniformly,  then  S  =  0.3  S  . 

u  c 


< 
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TABLE  5.1 

Containment  Scattering  Source  Strength 

Gamma  Energy  Total  Containment 

E  Source  Strength 

,  S 

(Mev)  c 

(Photons /sec) 


0.03 

1.59  x  1012 

0.04 

2.89  x  109 

0.05 

1.02  x  1012 

0.06 

9.27  x  10 10 

0.08 

5.89  x  1013 

0.10 

4.68  x  1010 

0.15 

3.18  x  1013 

0.20 

3.38  x  1013 

0.30 

3.42  x  1013 

0.40 

3.48  x  1013 

0.50 

2.51  x  1013 

0.60 

2.11  x  1013 

0.80 

5.10  x  1013 

1.0 

1.27  x  1013 

1.5 

1.75  x  1013 

2.0 

4.90  x  1013 

3.0 

5.59  x  1012 

4.0 

6.25  x  10^ 

The  dose  conversion  factor  requires  some  discussion.  Since 
the  dose  at  P  is  dependent  on  the  energy  of  the  scattered  photon, 

0^  in  this  case  is 

C  =  1.15  x  10“4  _  _  (E* ) 

D  a 


The  scattering  energy  E'  is  constantly  changin  ;  since  it  is  dependent 
on  the  angle  8.  Technically  then,  should  be  inside  the  summation 
and  calculated  each  time  a  new  8  (and  E')  is  determined.  This  comp¬ 
lication  can  be  simply  overcome  by  noting  that 


and 


where 


do 
_ s 

dft 


=  Klein-Nishina  differential  scattering  (energy) 
cross  section. 


which  is  given  by 


da 
_ s 

dfi 


(f)  (§,  +  f  -  sin2e)  cm2/steradian 


(5.10) 


The  collision  cross-section  calculates  the  fraction  of  the  flux 
which  is  scattered  into  8  while  the  scattering  cross-section  calculates 


the  fraction  of  incident  energy  scattered  into  6.  Therefore, 
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in  Eq .  5.9 


can  be  replaced  by 


da 

EdjT  ’ 


Since  E  is  a  constant  it  can  be  removed  from  the  integral.  The 
absorption  coefficient  y  (E')  can  be  conservatively  treated  by  replacing 

d 

it  with  y  (E) ,  which  is  defined  as  the  largest  value  of  y  (E’) 

3.  a 

within  the  scattered  energy  range  E  to  E'  (corresponding  to  180° 

min 

scattering).  Values  of  E'.  and  y  are  listed  in  Table  A. 6. 

mm  a 

Finally,  since  the  photons  must  pass  through  the  steel  shell 
credit  can  be  taken  for  attenuation  in  the  steel.  It  is  conservative 
to  use  the  minimum  (3/8  in)  thickness  b  (where  b  =  y  T„  is  the 

X  X  j  1  j  1 

steel  thickness  in  number  of  mean  free  paths.) 


The  two  hour  dose  at  point  P  is  then  given  by 

_b  TM>0  TT-^ 

1.15  X  10  4s  NEy  e  1  f  da 

Dose  = - - - -  /  A*  J  (6 


(6  =  ip  +  4> )  d<J>  rads  (5.11) 


The  above  equation  was  evaluated  for  each  energy  E  using  a 
simple  numerical  program  with  sufficient  intervals  (10)  to  yield 
no  change  in  the  third  significant  figure.  The  geometry  parameters 
used  are  shown  in  Table  A. 11  and  the  resulting  total  dose  at  each 
distance  due  to  all  energies  is  listed  in  Table  5.2. 


5.3.2  Lower  Containment 

Since  the  scattering  angle  changes  significantly  as  the  source 
moves  down  into  the  shadow  shield  cylinder,  the  lower  containment  will 


be  treated  as  three  equal  source  points  located  along  the  central 
axis  of  the  cylinder  at  a  height  equal  to  the  center  of  the  correspond 
ing  slice  of  the  containment.  These  points  will  be  designated 
1,  2,  and  3  as  shown  in  Fig.  5.3. 

In  this  case,  the  lateral  sweep  of  dV  in  the  cross-wise  direction 
is  limited  by  the  shadow  shield  to  some  value  less  than  it.  This 
azimuthal  angle  tu  is  dependent  on  the  depth  of  the  source  within 
the  shadow  shield  and  the  value  of  ip. 

Referring  to  Fig.  5.4a,  it  can  be  seen  that 

“  * 2  cos"1(f') 

Fig.  5.4b  provides  the  relationship 

sin  OP  -¥>)  =  Ip- 


Substituting,  the  result  is 


uj  -  2  cos 


-1 


R'  sin  (\p  -  f ) 


radians 


(5.12) 


For  a  given  point  along  the  central  axis  h' ,  R',  and  ip  are  constants. 
Therefore  a>  is  only  a  function  of  lji  and  can  be  placed  within  the 
scattering  integral  (in  this  case  within  the  summation)  in  place  of 
the  original  azimuthal  angle  tt. 


The  strength  of  each  point  is 


therefore  Eq.  5.11  becomes 


for  lower  containment  scattering  points 
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TT-<t>o  TT-lp 

1.15  X  10'4  S  NEy  e  1  V”1  f  do 

Dose  =  - 12ttx - - - -  /  Clp ) Ai pj  (0  =  ip  +  <p)d<p  rads 

4>0 

(5.13) 

where  S„  =  0.7S 
i  c 

and  the  values  x,  i|>Q ,  <J) Q ,  h',  R'  and  'f  are  dependent  on  the  source 

and  detector  geometry. 

Eq.  5.13  was  evaluated  for  each  source  point  and  each  energy 
using  the  parameters  in  Table  A. 11  and  the  resulting  doses  are  given 
in  Table  5.2. 

5.3.3  Air  Scattering  Results  and  Comparison  with  the  SAR 

The  total  air  scattering  doses  by  energy  level  are  given  in 
Table  5.3.  Note  that  the  more  isotropic  lower  energy  gammas  produce 
a  higher  dose  at  8  m  while  the  forward  scattering  high  energy  gammas 
yield  a  higher  dose  at  21  meters.  The  conservative  air  scattering 
model  used  in  this  section  -Stimates  the  total  two  hour  dose  as 
1.14  x  10  1  rads  at  8  meters  and  1.47  x  10  ^  rads  at  21  meters. 

The  total  21  meter  dose  is  larger  as  a  result  of  the  enhanced 
forward  peaking  due  to  the  smaller  scattering  angles  required  to 
reach  the  21  meter  point.  As  distance  continues  to  increase  the 
dose  will  peak  and  then  fall  off  as  the  x  term  begins  to  have  an 
effect.  An  estimate  of  the  peak  distance  and  dose  has  been  made 
by  tracking  the  behavior  of  the  two  highest  dose  energies,  E  =  0.8 
and  2.0  Mev,  for  each  source  point.  The  results,  based  on  the  ratio 


AD-R128  285 


UNCLASSIFIED 
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TABLE  5.2 

Air  Scattering  Doses  by  Source  Point 

Dose 

Source  Point  (Rads) 


8m 

21m 

s 

u 

6 . 89  x 

io-2 

9.98  x 

io~2 

.  si 

2.30  x 

Kf2 

2.46  x 

io'2 

S2 

1.34  x 

io"2 

1.39  x 

io"2 

s. 

8.32  x 

io-3 

8.55  x 

io'3 
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TABLE  5.3 

Total  Air  Scattering  Doses  by  Gamma  Energy 


Gamma  Energy 

Dose 

E 

(Mev) 

8m 

(Rads) 

21m 

-28 

0.03 

5.56 

X 

10 

5.16 

X 

10 

0.04 

4.00 

X 

IQ'17 

3.74 

X 

10 

0.05 

1.97 

X 

io"9 

1.85 

X 

10 

0.06 

2.04 

X 

10-8 

1.93 

X 

10 

0.08 

7.42 

X 

io'4 

7.08 

X 

10 

0.10 

2.30 

X 

10  “6 

2.23 

X 

10 

0.15 

5.63 

X 

io-3 

5.60 

X 

10 

0.20 

9.21 

X 

10“  3 

9.44 

X 

10 

0.30 

1.27 

X 

10' 2 

1.37 

X 

10 

0.40 

1.44 

X 

io'2 

1.70 

X 

10 

0.50 

1.07 

X 

io~2 

1.26 

X 

10 

0.60 

9.08 

X 

io"3 

1.11 

X 

10 

0.80 

2.15 

X 

io~2 

2.82 

X 

10 

1.0 

5.20 

X 

10"  3 

7.19 

X 

10 

1.5 

6.48 

X 

io-3 

1.01 

X 

10 

2.0 

1.61 

X 

io'2 

2.82 

X 

10 

3.0 

1.52 

X 

10"  3 

3.10 

X 

10 

4.0 

1.43 

X 

10~5 

3.31 

X 

10' 

Total 

1.14 

X 

io"1 

1.47 

X 

10' 

96 


of  these  components  to  the  total  dose,  are  shown  in  Fig.  5.9  and  indicate 
the  maximum  total  dose  is  0.1.5  x  10  ^  rads  at  a  distance  of  o.  27  meters. 

Only  single  scattering  events  have  been  considered.  This  is 
sufficient  as  it  has  been  shown  [34]  that  the  single  scattering 
model  derived  in  Sec.  5.2  yields  results  which  are  in  almost 
exact  agreement  with  Monte  Carlo  air  scattering  results  [35]. 

The  SAR  gives  air  scattering  values  of 0*1.55  rads  at  8  meters 
and  o,  8.3  x  10  rads  at  21  meters.  The  difference  is  due  to  the 
very  conservative  scattering  model  used  in  the  SAR,  which  does  not 
take  full  credit  for  the  effects  of  the  scattering  angles,  energy 
degradation  and  shadow  shield. 

5.4  Steel  Shell  Scattering 
5.4.1  Single  Scattering  Model 

The  air  scattering  model  could  be  directly  adopted  to  determine 
scattering  from  the  steel  shell  if  the  integration  could  be  made 
to  follow  the  surface  of  the  shell.  Since  this  is  not  easily  done 
the  following  approximation  will  be  made. 

If  the  entire  scattering  volume  shown  in  Fig.  5.5  were  composed 
of  steel,  then  the  dose  at  P  due  to  a  source  S  photons/sec  would  be 


‘  > 
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x 


(5.15) 


Substituting  in  Eq.  5.14,  the  dose  becomes 


Dose  * 


1.15  x  IQ-4  SNSIVsrEj;ae-b 

411^2 (*2  ”  ♦1)W2  -  *2) 


Alp  /  (6 ) dtp  rads 


(5.16) 


where  the  distances  r^  and  r^  are  shown  in  Fig.  5.6. 
The  electron  density  in  steel  is  determined  from 


where  for  steel  (iron) : 


3  23 

(26  elec/atom) (7 .83  g/cm  )(6.02  x  10  atoms/g-atom) 
55.85  g/g-atom 


nst  =  x  elec/cm^ 

The  volume  of  steel  in  the  dome  is  determined  from  Fig.  5.7  as 


VST  '  TST<0dLH  +  W 


9d  ’  2  Sl"_1  (2t> 

a 


therefore , 


VST  =  TST  2sln  1  (2R~)  (LH  +  Rdy)  Cm3 
„  d 


(5.17) 


where  =  Steel  thickness  =  9.5  x  10  cm  * 

ST 


£  =  Dome  segment  length  (cm) 


Although  the  dome  portion  of  the  shell  is  5/8  in  thick,  3/8  in 
was  used  here  for  the  entire  shell.  The  resulting  smaller  steel 
volume  is  partially  offset  by  the  decreased  attenuation.  The 
error  introduced  bv  using  3/8  in  is  estimsted  to  be  less  than  10% 


I 


Figure  5.7 

Steel  Shell  Volume  Parameters 


R  =  Dome  radius  =  1.13  x  10  cm 
d 

3 

L  =  Dome  curvature  =  2.13  x  10  cm 
H  =  Dome  vertical  distance  (cm) 

2 

y  =  Side  panel  height  =  1.68  x  10  cm 

Using  Eq.  5.17,  the  steel  volumes  are 
Vst(8)  =  8.19  x  105  cm3 

vst(21)  =  2,91  x  1q6  Cffl3 


5.4.2  Single  Scattering  Results  and  Comparison  with  the  SAR 

The  exclusion  area  boundary  doses  were  computed  using  Eq.  5.16  for 
the  same  four  source  points  as  used  in  air  scattering.  The  geometry 
parameters  for  each  point  are  given  in  Table  A. 12.  As  in  air  scatter¬ 
ing,  doses  due  to  photons  of  E  <  0.08  Mev  are  negligible  and  are 
not  included  in  the  total  steel  scattering  results  listed  in  Tables 
5.4  and  5.5.  The  estimated  total  dose  is  1.60  x  10  ^  rads  at  8  meters 
and  3.66  x  10  rads  at  21  meters.  The  larger  amount  of  steel  visible 
at  21  meters  causes  a  higher  dose. 

The  SAR  predicts  much  smaller  steel  scattering  doses,  namely 

-2 

0  rads  at  8  meters  and  n,  6.3  x  10  rads  at  21  meters. 

The  difference  can  be  explained  as  follows.  The  steel  scattering 


model  used  in  the  SAR  is  very  similar  except  that  the  scattering  from 
steel  is  assumed  to  be  isotropic  (as  opposed  to  differential).  There¬ 
fore,  for  a  given  source  strength,  location,  and  energy,  the  SAR  method 
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will  predict  a  higher  dose.  The  difference  in  results  lies  in  the  source 
strength  and  location.  As  stated  in  Chapter  2,  the  source  strength 
in  this  study  is  about  1.5  times  that  of  the  SAR.  The  difference 
in  location  is  that  the  SAR  used  a  single  source  point,  located 
in  the  center  the  shadow  shield.  The  effect  of  this  is  illustrated 

by  the  following  dose  equations  for  energy  E  =  0.8  Mev  and  a  distance 
of  21  meters  in  which  the  source  strength  has  been  left  as  S. 

From  the  SAR,  the  dose  is 

Dose  =  1.81  x  10  15  S  rads. 

In  this  study,  the  doses  are: 

Upper  Containment  Dose  =  3.45  x  10  ^  S 

Lower  Containment, 

Point  1  Dose  =  9.15  x  10  ^  S 

—16 

Point  2  Dose  =  3.60  x  10  S 

Point  3  Dose  =  1.92  x  10  ^  S 

The  SAR  source  point  corresponds  to  point  2  in  the  lower  con¬ 
tainment.  As  seen,  the  SAR  method  will  yield  a  higher  dose.  But 
as  the  source  point  moves  up  the  containment,  the  dose  increases. 

This  is  to  be  expected,  due  to  the  increase  in  forward  scattering  and 
the  decrease  in  scattering  distance.  The  dose  from  the  upper  containment 
is  almost  10  times  that  of  point  2.  This  sensitivity  to  location  combined 
with  the  higher  source  strength  results  in  higher  doses  than  the  SAR. 
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TABLE  5.5 

Total  Steel  Scattering  Doses  by  Gamma  Energy 


Gamma  Energy 
E 

(Mev) 


Dose 

(Rads) 

8m  21m 


0.08 

4.13 

X 

10-4 

6.82 

X 

1 

c 

rH 

0.10 

1.39 

X 

io-6 

2.28 

X 

IO-6 

0.15 

4.01 

X 

io"3 

6.64 

X 

io-3 

0.20 

7.45 

X 

10~3 

1.25 

X 

IO-2 

0.30 

1.25 

X 

io"  2 

2.19 

X 

io"2 

0.40 

1.63 

X 

io"  2 

2.97 

X 

io"2 

0.50 

1.35 

X 

io"2 

2.57 

X 

io"2 

0.60 

1.24 

X 

io-2 

2.44 

X 

io"2 

0.80 

3.35 

X 

io'2 

7.13 

X 

io-2 

1.0 

8.78 

X 

10~  3 

1.99 

X 

io"2 

1.5 

1.27 

X 

10"  2 

3.40 

X 

io"2 

2.0 

3.50 

X 

io"2 

1.06 

X 

io"1 

3.0 

3.69 

X 

io"3 

1.37 

X 

io-2 

4.0 

3.62 

X 

io"5 

1.42 

X 

io-4 

Total 

1.60 

X 

io" 1 

3.66 

X 

io-1 

The  importance  of  double  scattering  in  steel  can  be  evaluated 


as  follows. 

As  seen  in  Fig.  5.8  the  photons  which  impinge  on  the  dome  be¬ 
tween  points  A  and  B  can  only  arrive  at  P  due  to  a  second  scatter 
between  points  B  and  C.  To  be  conservative  it  is  assumed  that  all 
the  second  scatters  occur  at  point  C  as  this  increases  the  dose 
at  P  due  to  preferential  forward  scattering.  The  dose  at  point  C 
due  to  scatters  between  A  and  B  can  be  computed  using  the  single 
scattering  equation  as  before.  The  dose  which  scatters  from  C  to 
P  is  determined  by  multiplying  by  the  factor 


N  V'  do 

2  dfl  (6  ' 


where 


V'  =  Volume  of  steel  between  B  and  C. 


0'  =  Second  scattering  angle  =  <j>_  -  <|> 

Z  n 


the  resulting  two  scatter  dose  is 
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The  large  brackets  are  used  to  indicate  that  during  the  numer 


da 


ical  integration  for  each  value  of  <f>  both  _ s 


and  -777—  (0')  are  evaluated  and  their  product  determined. 
dW 

The  effect  of  double  scattering  has  been  estimated  by  evaluating 
Eq.  5.18  for  the  three  energies  (E  =  0.4,  0.8,  and  2.0  Mev)  which 
contribute  the  most  to  the  total  steel  dose  and  the  same  four  source 
points.  The  results  indicate  the  total  steel  scattering  dose  results 
should  be  increased  by  a  factor  of  1.20  at  8  meters  and  by  1.02 
at  21  meters.  The  adjusted  doses  are  1.92  x  10  rads  at  8  meters 
and  3.73  x  10  ^  rads  at  21  meters. 

The  behavior  of  the  steel  scattering  dose  past  21  meters  has 
been  investigated  by  tracking  the  behavior  of  the  dose  due  to  the 
same  three  energies  and  assuming  that  the  contribution  of  double 
scattering  past  21  meters  is  negligible. 

The  results,  shown  in  Fig.  5.9,  indicate  that  the  maximum 
steel  scattering  dose  is  ^  3.80  x  10  1  rads  at  a  distance  of  n,  22 
meters. 

5 . 5  Summary 

The  differential  scattering  model  developed  in  this  section 
and  applied  to  single  scattering  in  air  and  single  and  double  scattering 
in  steel  predicts  an  estimated  2  hour  dose  at  the  exclusion  area 
boundaries  of  3.06  x  10  ^  rads  at  the  back  fence  and  5.20  x  10 


rads  at  the  Albany  St.  fence.  The  estimated  behavior  of  the  total 
scattering  dose  is  shown  in  Fig.  5.9. 


Containment  (M) 
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CHAPTER  6 

RADIATION  PENETRATION  THROUGH  THE  TRUCK  LOCK 


6.1  General 

The  Integrity  of  the  shadow  shield  is  broken  by  a  number  of 
penetrations  for  piping,  cables  and  so  on,  the  majority  of  which  are 
small.  The  larger  penetrations;  the  truck  and  personnel  locks, 
the  neutron  window  and  the  hot  plug  storage  area;  are  shielded  by 
additional  concrete  walls  or  earth  [1].  The  largest  penetration 
is  the  truck  lock,  and  since  it  is  shielded  by  concrete  only  on  the 
sides,  the  possibility  exists  that  gamma  radiation  could  stream 
out  of  the  lock  and  reach  the  boundary  fence  by  either  reflection 
off  the  opposite  concrete  wall  or  scattering  through  the  steel  door 
(see  Fig.  6.1).  The  exclusion  area  dose  due  to  these  effects  will 
be  determined  in  the  following  sections. 

6.2  Truck  Lock  Description 

As  shown  in  Figs.  6.1  and  6.2,  the  truck  lock  is  a  rectangular 
steel  tube  8  meters  long  closed  at  both  ends  by  pneumatically  sealed 
doors  3  meters  wide  and  4  meters  high.  Each  door  consists  of  a 
steel  framework  covered  on  both  sides  by  0.6  cm  thick  steel  plates. 
The  inner  door  is  aligned  just  inside  the  containment  shadow  shield. 
The  two  sides  of  the  lock  are  shielded  by  concrete  walls  0.5  meters 


Figure  6.1  Truck  Lock  Radiation  Release  Pathways 


L2 


thick  while  the  front  and  top  are  unshielded  [1]. 


6.3  Truck  Lock  Source  Term 


To  facilitate  subsequent  calculations,  the  radiation  reaching 
the  lock  will  be  treated  as  a  point  source  located  at  the  center 
of  the  inner  surface  of  the  inner  door. 

The  amount  of  radiation  reaching  the  door  can  be  estimated 
by  using  the  same  four  containment  source  points  as  in  Chapter  5. 
If  each  point  emits  isotropically  then  the  number  of  photons/sec 
at  the  door  is  (see  Fig.  6.2) 
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where 


S  =  Truck  lock  source  strength  (photons/sec) 


T 


=  Area  of  the  truck  lock  door  (cm  ) 


As  before, 
and 


S 

u 


0.3  S 

c 


S„  0.7  S 

x,  _  _ c 

3  3 


where  S  =  S  +  S.  =  Total  containment  source  strength  (photons/sec) 

C  U  X  or/ 

The  source  distances  are: 


r  =  1.36  x 

10"^  cm 

r  =  1.07 

x  10^  cm 

u 

r,  =  1.15  x 

10^  cm 

2 

r_  =  1.09 

3 

x  10  cm 

I 
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The  truck  lock  area  is 

A^,  =  1.2  x  10^  cm^ 

Substituting,  the  final  truck  lock  source  term  is 

-3 

=  7.12  x  10  Sc  photons/sec  (6.2) 

where  the  values  of  are  given  in  Table  5.1. 

Finally,  it  is  conservatively  assumed  that  all  photons  reaching 
the  door  travel  through  it,  making  the  truck  lock  point  source  iso¬ 
tropic  in  the  forward  direction. 

One  could  argue  that  the  treatment  of  the  containment  source 
as  four  points  along  the  center  axis  does  not  give  full  weight  to 
those  fission  products  suspended  just  inside  the  door.  However,  it 
is  felt  that  this  treatment  is  still  conservative  since  the  entire 
containment  source  term  is  used,  while  in  reality  a  number  of  photons 
are  blocked  from  reaching  the  door  by  the  reactor  vessel.  In  addition, 
no  energy  degredation  is  considered,  although  the  majority  of  the 
gammas  would  have  to  make  a  change  of  direction  in  order  to  reach 
the  outer  truck  lock  door. 

6.4  Concrete  Scattered  Dose 

6.4.1  Unattenuated  Flux 

The  flux  on  the  concrete  wall  opposite  the  truck  lock  is  determined 
from  the  point  source  relationship  [30] 


A 


Th e  above  equation  neglects  build  up  and  attenuation  in  the  air 


The  2  in  the  denominator  (instead  of  4)  reflects  the  conservative 
assumption  that  the  source  is  isotropic  in  the  forward  direction. 
The  dose  at  the  wall  is  then 


1.15x10  BS^  _^T 

Dose  =  - 2 -  e 

2itx 


(6.4) 


The  minimum  thickness  of  steel  through  both  doors  is  2.54  cm 

3 

(1  in)  and  the  distance  to  the  wall  is  2.21  x  10  cm.  From  Eq.  6.4 
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TABLE  6.1 

Direct 

Dose  at  the  Concrete  Wall 

Gamma  Energy 

Build-up 

E 

Factor  EuT 

6.63 

1.54 

X 

10 

3.58 

3.44 

X 

10 

2.72 

1.24 

X 

10 

2.09 

3.19 

X 

10 

1.80 

5.59 

X 

10 

1.62 

5.22 

X 

10 

1.48 

5.59 

X 

10 

1.30 

1.89 

X 

10 

1.17 

5.98 

X 

10 

0.947 

1.23 

X 

10 

0.828 

4.39 

X 

10 

0.701 

6.48 

X 

10 

0.645 

8.76 

X 

10' 

Total 

1.09 

X 

10' 

Figure  6.3 
Albedo  Geometry 
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where  =  Energy  and  material  dependent  constants. 

The  total  dose  is  determined  by  integrating  over  the  area 
of  reflection.  The  area  of  the  wall,  A^,  in  view  of  the  source  is 
limited  by  the  truck  lock  side  walls  as  shown  in  Fig.  6.4.  Integra¬ 
tion  of  Eq.  6.5  is  difficult  since  the  angles  and  distances  are 
different  for  each  point  on  the  wall.  As  an  approximation,  the 
albedo  at  the  center  of  the  scattering  area  can  be  determined  and 
used  for  the  entire  area.  By  trial  and  error  the  point  along  the 
northern  boundary  was  found  which  receives  the  maximum  dose.  This 
is  point  in  Fig.  6.4  and  the  scattering  parameters  are: 

80  =  0.384  radians 

0  =  1.117  radians 

r 

9  =  1.641  radians 
r  =  2.17  x  10^  cm 

5  2 

The  area  of  A^  of  the  wall  in  view  of  the  source  is  4.18  x  10  cm  . 
Substituting  in  Eq.  6.5,  the  dose  is 
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D  =  8.23  x  10'2  DA  rads  (6.7) 

J  A. 

Using  the  values  of  and  for  concrete  given  in  [27]  a  graph 

of  A v  vs.  E  was  constructed.  The  resulting  values  of  A  and  the 
JX  JX 

two-hour  doses  are  shown  in  Table  6.2.  The  estimated  maximum  dose 
at  the  exclusion  area  boundary  due  to  backscattering  from  the 
concrete  wall  is  8.25  x  10  ^  rads. 


6.5  Steel  Door  Scattered  Dose 


The  second  path  along  which  gamma  radiation  can  reach  the 
exclusion  area  through  the  truck  lock  is  by  scattering  in  the.  outer 
steel  door.  The  steel  scattering  equation  developed  in  Chapter  5 
can  be  applied  directly  to  the  geometry  shown  in  Fig.  6.5.  The 
dose  at  is  then 
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1.15  x  10  S_N  V  Ep  e  ^ 

D  =  - — - T-S-T_ST__a -  )  ^ 
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do 
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(6.8) 


where 


-3 

S„  =  7.12  x  10  S  photons/sec 
i  c 

24  3 

N„„  =  2.19  x  10  elec/cm 
ST 


The  location  along  the  boundary  which  receives  the  maximum  dose  is 
point  P 2  with  the  parameters: 
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TABLE  6.2 

Concrete  Albedo  Dose 


Gamma  Energy 
E 

(Mev) 


Exposure  Albedo 
A ... 


Dose 

(Rads) 


0.10 

3.9 

X 

10-2 

A. 9  A 

X 

io"11 

0.15 

3.05 

X 

io"2 

8.64 

X 

10" 7 

0.20 

2.55 

X 

io"2 

2.60 

X 

io"6 

0.30 

2.0 

X 

10“  2 

5.26 

X 

10"6 

O 

o 

1.65 

X 

io"2 

7.61 

X 

io"6 

o 

m 

o 

1.A5 

X 

io"2 

6.22 

X 

10“6 

0.60 

1.30 

X 

io-2 

5.98 

X 

io"6 

0.80 

1.08 

X 

IO-2 

1.68 

X 

io-5 

1.0 

9. A 

X 

io'3 

4.63 

X 

io"6 

1.5 

7. A 

X 

10~3 

7.50 

X 

io"6 

2.0 

6.2 

X 

io"3 

2.24 

X 

io-5 

3.0 

A. 8 

X 

io"3 

2.56 

X 

io"6 

A  .0 

4.1 

X 

10"  3 

2.96 

X 

io"8 

Total 

8.25 

X 

io"5 

Figure  6.5 

Steel  Door  Scattering  Geometry 
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r  =  9.14  x  102  cm 
=  2.82  x  102  cm 

=  0.140  radians  =  0.436  radians 

<j>2  =  0.209  radians  i}>  =  0.768  radians 

The  volume  of  steel  in  one  door  was  determined  from  the  dimensions 

and  material  list  in  the  original  reactor  housing  plans  [36]  to 
5  3 

be  2.01  x  10  cm  .  As  in  the  direct  dose,  an  attenuation  of  1  inch 
of  steel  was  used,  attenuation  and  build  up  in  the  air  have  been 
neglected,  and  the  factor  of  2  in  the  denominator  reflects  the 
forward  direction  source. 

As  shown  in  Table  6.3,  the  maximum  estimated  dose  due  to  steel 
■“A 

scattering  is  3.55  x  10  rads. 

6.6  Air  Scattered  Dose 

Using  the  techniques  of  Chapter  5  the  effect  of  air  scattering 
from  the  truck  lock  to  the  exclusion  area  was  investigated  and  found 
to  be  less  than  1  x  10  2  rads.  This  is  primarily  due  to  the  limited 
scattering  angles  as  a  result  of  the  concrete  side  walls  along 
the  truck  lock. 

6.7  Conclusions 

The  concrete  and  steel  scattering  results  show  that  the  est¬ 
imated  radiation  release  from  the  reactor  truck  lock  during  the 
two-hour  period  does  not  add  significantly  to  the  total  dose  at  the 
exclusion  area  boundary  from  the  other  modes  of  release. 


CHAPTER  7 


SUMMARY,  RESULTS  AND  CONCLUSIONS 


7.1  Summary 

7.1.1  Introduction 


The  purpose  of  this  thesis  was  to  determine  a  conservative 
estimate  of  the  maximum  amount  of  radiation  an  individual  located 
at  the  exclusion  area  boundary  of  the  MIT  reactor  would  receive 
during  the  first  two  hours  of  the  reactor's  Design  Basis  Accident 
(DBA) . 

The  DBA  for  the  MITR  is  postulated  to  be  a  coolant  flow  blockage 
in  one  of  the  flat  plate  fuel  elements  resulting  in  complete  melting 
of  four  U-Al^  fuel  plates.  The  radiation  release  from  this  event 
has  been  previously  examined  in  the  MIT  Reactor  Safety  Analysis 
Report  (SAR) ,  published  in  1970. 

The  present  study  offers  a  significant  improvement  over  the 
SAR  in  the  following  areas: 

1)  The  fission  product  inventory  has  been  expanded  to  include 
other  important  isotopes,  such  as  those  of  Tellurium. 

2)  The  release  fractions  have  been  updated  to  more  realistic 
values  based  on  the  expected  physical  processes  specific 
to  the  MITR. 

3)  The  leakage  release  has  been  based  on  historical  meteoro¬ 
logical  data  and  includes  gamma  and  beta  whole-body  doses. 


4)  The  scattering  model  of  gamma  radiation,  particularly  from 
air,  has  been  greatly  improved  and  is  no  longer  grossly 
over-conservative. 

5)  The  effect  of  penetrations  have  been  examined  for  the  first 
time . 

7.1.2  Development  of  the  Containment  Source  Term 

The  fission  product  isotopes  of  interest  include  the  noble 
gases,  iodine  and  any  less  volatile  element  which  can  be  expected 
to  yield  a  measurable  whole-body  or  thyroid  dose.  The  fission  product 
inventory  in  the  fuel  at  the  time  of  the  accident  is  assumed  to  be 
equal  to  the  maximum  equilibrium  value  after  an  essentially  indefinite 
period  of  irradiation  at  a  power  level  of  5  MW.  The  saturation 
activities  of  the  fission  products  can  be  determined  using  analytical 
production  and  decay  relationships  or  numerical  results  such  as  those 
by  Blomeke  and  Todd  [8].  Once  the  total  core  fission  product  inventory 
is  known,  the  fraction  released  along  each  portion  of  the  path  to 
the  containment  must  be  determined. 

Assuming  a  radial  peaking  factor  of  1.45  the  fraction  of  the 
core  fission  product  inventory  contained  in  the  four  melted  fuel 
plates  is  1.5% 

Assuming  typical  End-of-Cycle  fuel  conditions,  studies  of 
U-Al  melting  [6,  12,  13,  14]  show  that  the  estimated  fractions  of 
the  fission  products  in  the  melted  fuel  released  to  the  primary 


coolant  are: 
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100  %  of  the  noble  gases  (Kr,  Xe) 

100  %  of  the  halogens  (I,  Br) 

70%  of  the  Tellurium 

30%  of  the  alkali  metals  (Cs,  Rb) 

1%  of  the  remaining  fission  products 

Release  from  the  coolant  water  is  dependent  on  the  degree  of 
plateout  within  the  primary  system  and  absorption  in  the  coolant . 

Recent  studies  [5,  16,  19]  indicate  that  iodine  will  be  retained  in 

the  coolant  in  the  form  of  Csl  to  a  much  greater  degree  than  previously 

expected,  with  retention  as  high  as  99.9%.  For  the  present, 
however,  the  following  conservative  (especially  with  respect  to 
iodine)  values  will  be  used: 

100%  of  the  noble  gases 
10%  of  all  other  isotopes. 

Natural  depletion  in  the  containment  atmosphere  is  dependent 
on  deposition  velocities.  Based  on  the  findings  in  WASH-1400  [15], 
the  fraction  of  the  released  fission  products  which  remain  airborne 
are  expected  to  be: 

100%  of  the  noble  gases 
30%  of  the  halogens 

90%  of  the  remaining  fission  products. 

The  resulting  overall  airborne  release  fractions  compare 
favorably  with  those  of  historical  accidents,  while  still  being 
conservative  with  respect  to  the  latest  source  term  revision  estimates 
U9]. 


7.1.3  Atmospheric  Release 


It  is  anticipated  that  fission  products  in  the  containment 
atmosphere  may  leak  out  of  the  containment  through  small  cracks 
and  crevices.  It  is  conservatively  assumed  that  this  occurs  at  the 
maximum  leakage  rate  of  1%  of  the  building  volume  per  psi  over¬ 
pressure  per  day  and  that  the  containment  pressure  is  at  the  max¬ 
imum  value  of  2  psi  over  atmospheric. 

Once  released,  the  concentration  of  fission  products  in  the  air 
will  be  decreased  due  to  atmospheric  diffusion,  which  is  modeled 
based  on  the  Gaussian  plume,  and  the  building  wake  effect  [22],  This 
wake  effect  reflects  the  tendency  for  the  released  particles,  to 
quickly  expand  to  fill  a  volume  of  air  on  the  order  of  the  contain¬ 
ment  volume.  The  relative  magnitude  of  these  two  effects  depends 
on  the  distance  downwind  from  the  containment.  In  the  range  of  the 
exclusion  area,  the  wake  effect  is  extremely  dominant. 

The  NRC  has  proposed  a  conservative  concentration  equation  which 
takes  only  partial  credit  for  the  wake  effect  [23].  Using  this 
equation  with  wind  data  from  Boston's  General  Logan  Airport  the 
estimated  whole-body  dose  due  to  gamma  and  beta  radiation  to  a  person 
located  at  the  minumum  exclusion  area  distance  of  8  meters  is  8.12  rads 
and  the  dose  to  the  thyroid  is  36.3  rads.  At  the  most  populated  edge 
of  the  exclusion  area  (a  distance  of  21  meters)  the  whole-body  dose 
is  1.23  rads  and  that  to  the  thyroid  is  5.49  rads. 

In  contrast,  the  estimated  dose  taking  into  account  the  full  wake 


effect  is  2.66  x  10  rads  whole-body  and  1.18  x  10  rads  thyroid 
at  any  point  on  the  exclusion  area  boundary.  The  actual  dose  is 
expected  to  be  somewhere  in  between  these  values.  Due  to  the 
conservative  assumptions  with  respect  to  the  building  pressure  and 
leakage  rate,  plus  the  wind  tunnel  effect  of  adjacent  buildings, 
the  lower  value  is  deemed  more  likely. 

7.1.4  Direct  Gamma  Dose 

Those  isotopes  which  do  not  leak  from  the  containment  will 

constitute  a  source  of  gamma  radiation.  A  circular  concrete  shadow 

shield  9.6  m  high  and  0.61  m  thick  is  incorporated  into  the  containment 

to  protect  against  this  hazard.  Above  the  shadow  shield  the  gammas 

are  attenuated  only  by  the  0.95  cm  thick  steel  containment  shell. 

From  any  point  on  the  exclusion  area  boundary  some  portion  of 

the  containment  volume  is  visible  above  the  top  of  the  shadow  shield. 

This  volume  can  be  modeled  as  a  spherical  volume  source  shielded 

by  a  steel  slab.  Using  available  point  kernel  integration  expressions 

[27]  with  the  Taylor  analytical  form  of  the  build  up  factor  the  dose 

-3 

at  8  meters  from  the  containment  is  estimated  to  be  3.49  x  10  rads 

-2 

while  that  at  a  distance  of  21  meters  is  2.71  x  10  rads. 

That  portion  of  the  containment  volume  located  behind  the  shadow 


shield  can  be  modeled  as  a  right  circular  cylinder  shielded  by  a 
two  layer  concrete/steel  slab.  Point  kernel  results  for  this  shape 
have  been  determined  by  numerical  methods  [27],  and  using  the  build  up 
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factor  for  concrete,  yield  an  estimated  two-hour  dose  of  4.33  x  10 

-2 

at  8  meters  and  2.05  x  10  at  21  meters. 

The  maximum  total  direct  dose  occurs  beyond  the  exclusion  area 

-2 

at  a  distance  of  'v  38  meters  with  a  value  of  ^  7.2  x  10  rads. 

7.1.5  Scattered  Gamma  Dose 

It  is  also  possible  for  gamma  radiation  to  reach  the  exclusion 
area  boundaries  by  escaping  over  the  top  of  the  shadow  shield  and 
scattering  back  down  from  interactions  with  the  air  or  the  steel 
containment  dome. 

Separate  scattering  models  for  air  and  steel  dome  scattering 
were  developed  based  on  the  Klein-Nishina  differential  scattering 
cross-section  which  take  into  account  the  scattering  angle  restrictions 
due  to  the  presence  of  the  shadow  shield.  The  volumetric  containment 
source  was  approximated  as  four  point  sources  located  along  the  central 
axis  of  the  containment.  Attenuation  in  the  steel  shell  and  photon 
energy  degradation  due  to  scattering  were  also  considered. 

The  resulting  expression  was  numerically  integrated  over  the 
range  of  scattering  angles  to  yield  an  estimated  air  scattering  dose 
of  1.14  x  10  1  rads  at  8  meters  and  1.47  x  10  ^  rads  at  21  meters. 
Including  double-scattering,  the  steel  dose  is  expected  to  be 
1.92  x  10  1  rads  at  8  meters  and  3.73  x  10  ^  rads  at  a  distance  of 
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21  meters.  The  maximum  total  scattering  dose  occurs  at  'v  23  meters 
and  is  believed  to  be  n-  5.3  x  10  ^rads. 

7.1.6  Radiation  Penetration  through  the  Truck  Lock 

The  effect  of  penetrations  in  the  shadow  shield  was  investigated 

by  calculating  the  dose  at  the  exclusion  area  due  to  travel  through 

the  largest  penetration,  the  truck  airlock.  There  are  three  pathways 

to  the  exclusion  area  boundary;  direct  radiation  backscattering  from 

the  opposite  concrete  wall,  scattered  radiation  from  the  outer  steel 

door,  and  scattered  radiation  from  the  air. 

The  direct  radiation  passing  through  the  doors  was  evaluated 

using  a  standard  attenuation  relation  with  a  steel  build  up  factor. 

Using  an  analytical  albedo  model  and  for  the  given  geometry  the 

-5 

maximum  boundary  dose  is  expected  to  be  8.25  x  10  rads. 

Steel  and  air  scattering  were  evaluated  using  the  previously 

developed  scattering  model  to  yield  a  maximum  steel  door  scattering 
-4 

dose  of  3.55  x  10  rads  and  an  air  scattering  dose  of  less  than 
1  x  10  rads. 

7.2  Results 


A  summary  of  the  estimated  doses  due  to  all  modes  of  radiation 
release  plus  the  resulting  total  release  are  shown  in  Table  7.1  and 
Fig.  7.1.  The  "exact"  (full  wake)  values  of  the  leakage  doses  have 
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been  used  since,  as  previously  mentioned,  they  are  believed  to  be 
more  realistic.  In  fact,  they  can  be  considered  conservative  with 
respect  to  a  realistic  building  pressure  and  leakage  points  distributed 
around  the  containment. 

The  conservative  assumptions  and  calculational  models  used  in 
this  investigation  give  an  estimated  maximum  dose  to  an  individual 
located  at  the  exclusion  area  boundary  during  the  first  two  hours 
of  the  MITR  design  basis  accident  of  5.95  x  10  ^  rads  whole-body  and 
1.18  x  10  1  rads  to  the  thyroid.  This  dose  is  the  integrated  total 
value  during  the  first  two  hours  following  release.  Fission  product 
decay  and  aerosol  deposition  will  greatly  reduce  the  dose  rate  after 
the  first  two  hours.  This  maximum  value  is  estimated  to  occur 
at  a  distance  of  21  through  23  meters  from  the  containment  wall, 
corresponding  to  the  exclusion  area  boundary  along  Albany  St. 

7.3  Conclusions 

Accordin  'o  10  CFR  100  [9]  the  maximum  calculated  dose  at 
any  point  along  the  exclusion  area  boundary  at  the  end  of  the  first 
two  hours  following  the  onset  of  fission  product  release  from  a  design 
basis  accident  must  be  less  than  25  rads  to  the  whole-body  and  300 
rads  to  the  thyroid.  Based  on  toe  above  results  it  can  be  concluded 
that  the  present  containment  design  and  exclusion  area  dimensions  for 
the  MIT  Reactor  are  well  within  the  regulatory  standards. 


TABLE  7.1 

Total  Dose  Results 


Figure  7.1  Total  Two  Hour  Whole-body  Dose  Results 


There  are  two  ways  in  which  the  results  of  this  study  could 
be  improved  upon  with  very  little  modification  of  the  calculational 
procedure.  They  are:  improved  release  fractions  and  additional 
scattering  source  points. 

During  the  period  in  which  this  report  was  prepared  a  number 
of  independent  studies  were  begun  by  various  organizations  with  the 
purpose  of  improving  on  the  various  fission  product  release  fractions 
reported  in  WASH-1400.  It  is  expected  that  the  results  of  these 
studies  will  show  that  actual  release  fractions  are  significantly 
lower  than  those  used  in  this  study. 

Although  the  scattering  model  developed  in  Chapter  5  has  resulted 
in  a  lower  scattering  dose  than  the  SAR,  it  is  still  the  largest 
component  of  the  total  dose.  As  shown  in  Sec.  5.4.2,  the  location  of 
the  gamma  point  source  within  the  containment  has  a  considerable 
effect  on  the  magnitude  of  the  dose  at  the  boundary.  For  simplicity,  only 
four  source  points  were  used  in  this  study,  with  all  four  located 
along  the  central  axis  of  the  containment.  Obviously,  greater 
accuracy  could  be  obtained  by  the  use  of  additional  source  points 
throughout  the  containment.  In  addition  this  would  allow  the 
shielding  and  scattering  angle  effects  of  the  reactor  to  be  included. 

One  area  which  was  not  addressed  in  this  thesis,  but  which  would 
be  of  interest  is  the  time-dependent  dose  rate  at  the  exclusion  area 
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NOMENCLATURE 


Parameter  of  the  direct  dose  volume 

Coefficient  of  the  Taylor  build  up  factor  expression 
Differential  exposure  albedo 

2 

Containment  surface  area  (cm  ) 

Atomic  weight  (g/g-atom) 

2 

Area  of  the  truck  lock  door  (cm  ) 

2 

Area  of  the  concrete  albedo  wall  (cm  ) 

2 

Containment  cross-sectional  area  (m  ) 


Shield  thickness  in  mean  free  paths 
Steel  dome  thickness  in  mean  free  paths 
Shadow  shield  thickness  in  mean  free  paths 
Build  up  factor 

3 

Breathing  rate  (m  /sec) 


Shape  factor 

2 

Flux  to  dose  conversion  factor  (rads  per  photon/cm  -s) 

3 

Initial  fission  product  concentration  in  containment  (Ci/m  ) 

3 

Fission  product  concentration  in  containment  at  time  t(Ci/m  ) 
Thyroid  dose  conversion  factor  (rads  per  Ci  inhaled) 

Sutton  lateral  diffusion  parameter 
Sutton  vertical  diffusion  parameter 


mi 
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NOMENCLATURE  (Cont.) 

Albedo  constant 

C2  Albedo  constant 

Cy  Photon  dose  conversion  factor  (rads  per  Ci-sec/m3) 


d 


d* 


D 


D 

o 


Parameter  of  the  direct  dose  volume  V 
Parameter  of  the  angle  to  (cm) 

Dose  (rads) 

Albedo  incident  dose  (rads) 

Thyroid  dose  (rads) 

Beta  dose  (rads) 

Beta  dose  rate  (rads/sec) 

Gamma  dose  (rads) 

Gamma  dose  rate  (rads/sec) 


e 

E 


E' 


Y 


Parameter  of  the  direct  dose  volume  V 
Photon  energy  (Mev) 

Scattered  photon  energy  (Mev) 

Special  exponential  integral  function 

Average  beta  energy  per  disintegration  (Mev/dis) 

Average  gamma  energy  per  disintegration  (Mev/dis) 


Fc  Fraction  of  fission  products  released  from  the  primary  system 

remaining  airborne  in  the  containment  atmosphere 

Ff  Fraction  of  fission  products  released  from  the  melted  fuel 

Fp  Fraction  of  fission  products  released  from  the  primary  system 


NOMENCLATURE  (Cont.) 


Overall  fraction  of  the  total  core  fission  product  inventory 
remaining  airborne  in  the  containment  atmosphere 

Fraction  of  the  total  saturated  core  fission  product 
inventory  present  in  the  melted  fuel. 

Cylindrical  source  geometry  attenuation  function 

Height  (cm) 

Parameter  of  the  angle  u>  (cm) 

Dimension  of  the  steel  dome  (cm) 

Isotope  i 

Parameter  of  the  G  function 

Dome  segment  length  (cm) 

Dome  curvature  (cm) 

Meander  correction  factor 

Sutton  parameter  associated  with  stability 

3 

Electron  density  in  air  (elec/cm  ) 

Number  of  nuclei  of  isotope  i 

235 

Initial  number  of  U  atoms  in  fuel 
Saturation  number  of  nuclei 

3 

Electron  density  in  steel  (elec/cm  ) 
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NOMENCLATURE  (Cont.) 
p  Parameter  of  the  G  function 

P  Reactor  power  (Mw) 

Q  Fission  product  release  rate  from  the  containment  (Ci/sec) 

Qc(t)  Activity  in  the  containment  at  time  t  (Ci) 

Q  Total  decay  emissions  in  the  containment  (dis) 

CT 

Qg  Fission  product  saturation  activity  (Ci) 

Q(t)  Fission  product  release  rate  at  time  t  (Ci/sec) 

Total  fission  product  release  (Ci) 

Q, j,(x)  Total  fission  product  release  at  distance  x  (Ci) 

r  Distance  (cm) 

r  Classical  radius  of  the  electron  (cm) 

e 

R  Radius  (cm) 

R'  Parameter  of  the  angle  oj  (cm) 

R,  Radius  of  the  steel  dome  (cm) 

d 

R^  Radius  associated  with  (cm) 

Radius  associated  with  (cm) 

s  Shadow  shield  penetration  distance  parameter  (m) 

S  Gamma  point  source  strength  (photons/sec) 

2 

SA  Surface  source  strength  (photons/cm  -sec) 

Sc  Containment  total  source  strength  (photons/sec) 

S^  Lower  containment  source  strength  (photons/sec) 
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NOMENCLATURE  (Cont.) 

Truck  lock  source  strength  (photons/sec) 

Upper  containment  source  strength  (photons/sec) 

3 

Volume  source  strength  (photons/cm  -sec) 

Time  (sec) 

Half-life 

Shield  thickness  (cm) 

Concrete  thickness  (cm) 

Steel  thickness  (cm) 

Windspeed 

Deposition  velocity  (cm/sec) 

3 

Containment  volume  (cm  ) 

3 

Volume  of  the  fuel  (cm  ) 

3 

Lower  containment  volume  (cm  ) 

3 

Steel  volume  (cm  ) 

3 

Volume  of  steel  subject  to  second  scattering  (cm  ) 

3 

Total  scattering  volume  (cm  ) 

3 

Upper  containment  volume  (cm  )' 

3 

Direct  dose  steel  penetration  volume  (cm  ) 

3 

Direct  dose  shadow  shield  penetration  volume  (cm  ) 

Distance  (cm) 


Sector  distance  (m) 
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NOMENCLATURE  ( Cont . ) 


Dome  side  panel  height  (cm) 
Fission  product  yield  (%) 

Self-absorption  distance  (cm) 
Atomic  number 


Coefficient  of  the  Taylor  build  up  factor  expression 
Coefficient  of  the  Taylor  build  up  factor  expression 

Parameter  used  in  determination  of  the  difference  between 
two  values  of  the  E^function. 

Compton  scattering  angle  (radians) 

Steel  second  scattering  angle  (radians) 

Dome  segment  angle  (radians) 

Albedo  incident  angle  (radians) 

Albedo  angle  of  reflection  (radians) 

Steel  penetration  disk  source  angle  (radians) 

Radiological  decay  constant  (sec  ^) 

Containment  leakage  rate  (V/sec) 

Containment  removal  rate  (sec  ^) 

Linear  attenuation  coefficient  (cm  ^) 
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NOMENCLATURE  (Cont.) 


True  energy  absorption  coefficient  in  air  (cm  /g) 


Maximum  air  absorption  coefficient  in  the  range  E  to 
(cm  /g) 


E' 

min 


Linear  attenuation  coefficient  for  concrete  (cm  ) 

Linear  attenuation  coefficient  for  the  source  medium  (cm  }) 
Linear  attenuation  coefficient  for  steel  (cm 


Density  (g/cm  ) 


Macroscopic  fission  cross-section  (cm  ) 


cr^  Microscopic  fission  cross-section  (cm  ) 

Oy  Lateral  Pasquill  plume  dispersion  coefficient  (m) 

a  Vertical  Pasquill  plume  dispersion  coefficient  (m) 

2 


do 

dQ 


do, 

dft 


Klein-Nishina  differential  collision  cross-section 
(cm^/steradian) 


Klein-Nishina  differential  scattering  cross-section 
(cm^/steradian) 


Scattering  angle  with  respect  to  the  detector  (radians) 

4>0  Initial  value  of  $  (air  scattering)  (radians) 

2 

4)^,  Thermal  neutron  flux  (neutrons/cm  -sec) 

4>^  Lower  limit  of  4)  (steel  scattering)  (radians) 

4^  Upper  limit  of  41  (steel  scattering)  (radians) 

2 

4)^  Photon  flux  (photons/cm  -sec) 


Parameter  of  the  angle  to  (radians) 

3 

Fission  product  concentration  in  air  (Ci/m  ) 

3 

Relative  concentration  (sec/m  ) 

Scattering  angle  with  respect  to  the  source  (radians) 
Initial  value  of  \p  (air  scattering)  (radians) 

Lower  limit  of  ip  (steel  scattering)  (radians) 

Upper  limit  of  ip  (steel  scattering)  (radians) 

Azimuthal  scattering  angle  (radians) 
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TABLE  A.  3 

Relative  Frequency  of  Stability  Conditions, 
Blue  Hills,  MA,  1888-1889  [2] 


Stability 

Condition 

Pasquill 

Category 

Relative 

Frequency 

Uns  table 

A,B ,  C 

30.5 

Neutral 

D 

24.5 

Stable 

E 

39.5 

Stable 

F 

5.0 

Stable 

G 

0.5 

TABLE  A. 4 
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TABLE  A. 5 

Gamma  Emission  Energies  by  Isotope 


Isotope 

Total  Decay 
Emissions 

T 

(Dis) 

Photon  Energy  (Mev) 
and  Distribution  (%)  [38,39] 

Kr 

85m 

2.22  x 

io17 

0.15  (78),  0.3  (14) 

87 

2.87  x 

10 17 

0.4  (50),  0.8  (8),  3.0  (14) 

88 

5.03  x 

io17 

0.03  (2),  0.15  (7),  0.2  (35),  0.4  (5), 

0.8  (23),  1.5  (14),  2.0  (53) 

Xe 

131m 

5.18  x 

io15 

0.15  (2) 

133m 

2.73  x 

io16 

0.2  (10) 

133 

1.12  x 

io18 

0.08  (37) 

135m 

5.80  x 

10 16 

0.5  (81) 

135 

1.53  x 

io17 

0.3  (91),  0.6  (3) 

138 

1.93  x 

io17 

0.15  (10),  0.3  (30),  0.4  (12),  0.5  (3) 

2.0  (37) 

I 

131 

4.97  x 

io16 

0.08  (3),  0.3  (5),  0.4  (82),  0.6  (7), 
0.8  (2) 

132 

5.76  x 

io16 

0.5  (20),  0.6  (99),  0.8  (85),  1.0  (22) 
1.5  (8),  2.0  (2) 

133 

1.09  x 

io17 

0.5  (86) 

134 

6.59  x 

10 16 

0.15  (3),  0.4  (8),  0.5  (8),  0.6  (18), 
0.8  (160),  1.0  (11),  1.5  (9),  2.0  (5) 

135 

9.19  x 

10 16 

0.4  (6),  0.8  (8),  1.0  (38),  1.5  (46), 

2.0  (10) 

Br 

83 

6 . 31  x 

10 15 

0.5  (1.4) 

84 

6 . 43  x 

10 15 

0.6  (1),  0.8  (48),  1.0  (8),  2.0  (25), 

4.0  (7) 

Cs 

134 

3.43  x 

io16 

0.5  (1),  0.6  (121),  0.8  (95),  1.0  (3), 

1.5  (3) 

» 


9 
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TABLE  A. 5  (Cont.) 

Isotope  Total  Decay  Photon  Energy  and  Distribution 

Emissions 


Cs 

136 

4.95 

X 

io15 

0.06 

(11) 

,  0.08  (6), 

0. 

15  (36) 

9 

0.3  (71), 

1016 

0.8 

(100) 

,  1.0  (82), 

1. 

5  (20) 

137 

2.77 

X 

0.6 

(85) 

Rb 

86 

7.32 

X 

10 15 

1.0 

(9) 

Te 

127m 

6.76 

X 

10 14 

0.06 

(1) 

127 

6.14 

X 

10 11 

0.4 

(1) 

129m 

4.10 

X 

10 15 

0.6 

(3) 

129 

7.17 

X 

io15 

0.03 

(17) 

,  0.5  (7), 

1.0 

(1) 

131m 

4.25 

X 

10 15 

0.08 

(2). 

0.10  (5), 

0.2 

(16), 

0. 

3  (9), 

io16 

0.8 

(91), 

1.0  (24), 

1.5 

(3),  2 

.0 

(3) 

131 

1.01 

X 

0.15 

(68) 

,  0.5  (21), 

0. 

6  (4), 

1. 

0  (13) 

132 

5.26 

X 

io16 

0.05 

(14) 

,  0.2  (88) 

133m 

3.09 

X 

io16 

0.3 

(ID, 

0.4  (1),  0 

.6 

(23),  0 

.8 

(8), 

io16 

1.0 

(89) 

134 

3.64 

X 

0.08 

(21) 

,  0.2  (48), 

0. 

3  (21), 

0 

.4  (19), 

0.5 

(35), 

0.8  (45) 

Sr 

91 

9.49 

X 

10 14 

0.6 

(15), 

0.8  (27), 

1.0 

(33), 

1. 

5  (5) 

Ba 

140 

1.08 

X 

10 15 

0.03 

(U) 

,  0.15  (6), 

0. 

3  (6), 

0. 

4  (5), 

0.5 

(34) 

Ru 

103 

4.99 

X 

io14 

0.5 

(88), 

0.6  (6) 

105 

1.34 

X 

10 14 

0.3 

(17), 

0.4  (6),  0 

.5 

(20),  0 

.6 

(16), 

io13 

0.8 

(48) 

106 

6.55 

X 

No  Gamma  : 

Decay 

RL 

105 

1.52 

X 

io14 

0.3 

(24) 

Tc 

99m 

9.24 

X 

io13 

0.15 

(90) 

Mo 

99 

1.04 

X 

1  ol3 

0.04 

(2), 

0.2  (7),  0 

.4 

(1),  o. 

8 

(16) 
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Isotope  Total  Decay 
Emissions 


TABLE  A. 5  (Cont.) 

Photon  Energy  and  Distribution 


Sb  127 


4.28  x  10' 


1.48  x  10 J 


0.06  (1),  0.3  (3),  0.4  (9),  0.5  (29) 
0.6  (45),  0.8  (17) 

0.4  (5),  0.5  (21),  0.6  (12),  0.8  (58), 
1.0  (46) 


Nd  147 


La  140 


Ce  141 


Zr  95 


4.46  x  10 


1.07  x  10 


1.03  x  10 
1.05  x  10] 


1.05  x  10J 


1.11  x  10 
1.03  x  10] 


0.10  (28),  0.3  (3),  0.4  (4),  0.5  (13) 

0.3  (20),  0.5  (40),  0.8  (19),  1.0  (10) 
1.5  (96),  3.0  (3) 

0.15  (48) 

0.06  (11),  0.3  (46),  0.5  (3),  0.6  (7), 

0.8  (10),  1.0  (1) 

0.08  (2),  0.15  (11) 

0.8  (98) 

0.6  (92) 


Nb  95 


1.11  x  10J 


0.6  (100) 
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TABLE  A.  7 

Shield  Thicknesses  in  Mean  Free  Paths 


Gamma  Energy 

E 

(Mev) 

Concrete 

UT  (b) 

Steel  (b^) 

Total  (b 

0.03 

160.4 

56.9 

217.3 

0.04 

79.9 

25.0 

104.9 

0.05 

51.7 

13.3 

65.0 

0.06 

39.2 

8.17 

47.4 

0.08 

28.7 

3.98 

32.7 

0.10 

24.5 

2.48 

27.0 

0.15 

20.1 

1.34 

21.4 

0.20 

17.9 

1.02 

18.9 

0.30 

15.3 

0.780 

16.1 

0.40 

13.7 

0.672 

14.4 

0.50 

12.5 

0.604 

13.1 

0.60 

11.6 

0.555 

12.2 

0.80 

10.1 

0.485 

10.6 

1.0 

9.15 

0.437 

9.57 

1.5 

7.50 

0.354 

7.85 

2.0 

6.41 

0.310 

6.72 

3.0 

5.23 

0.262 

5.49 

4.0 

4.58 

0.241 

4.82 

TABLE  A. 8 

Point  Isotropic  Source  Exposure  Build- 
Factors  for  Iron  (Steel)  [40,41] 
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APPENDIX  B 

CONTAINMENT  DIMENSIONS  AND  PARAMETERS 


A  cross-sectional  view  of  the  containment  with  its  dimensions  is 
shown  in  Fig.  B.l.  The  effect  of  smaller  equipment  located  within  the 
containment  on  the  volume  and  surface  area  will  be  neglected.  Only  the 
reactor  vessel  and  its  shielding  are  taken  into  account. 

B.l  Containment  Volume 


For  convienience  the  containment  can  be  considered  as  the  three 
sections  shown  in  Fig.  B.2.  Only  the  volume  above  grade  is  considered. 
This  is  conservative  as  it  results  in  a  higher  fission  product 
concentration  and  increased  source  strength.  During  an  actual  release 
some  fission  products  will  travel  into  the  basement. 

For  the  dome,  the  volume  is 

V-  ?  *  V3rl  +  h4> 

For  the  center  slice, 

V  =  tt  r^  h3 

and  for  the  lower  section, 

2  2 

V  =  tt  h2  -  tt  r^  h^ 

Substituting  values  from  Fig.  B.l,  the  total  volume  is 

V  =  1.67  x  105  ft3  (  4.73  x  103  m3  ) 


IMSUtM  \OH 


Figure  B.l  Cross  Section  of  Reactor  Building  [1] 
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B.2  Containment  Surface  Area 

For  the  dome  the  area  is 

A  =  2  it  h.  L 
4 

where  L  =  Radius  of  dome  curvature  =  70  ft. 

The  area  of  the  sides  is  equal  to 
A  =  2  it  r1(  h2  +  h3  ) 

The  area  due  to  the  floor,  reactor  and  shadow  shield  top  is 
2 

A  =  ir  r1  +  2  tt  r3  h1 

with  the  resulting  total  inner  surface  area  of 
Ag  =  1.91  x  10  4  ft2  (  1.78  x  103  m3  ) 

B.3  Containment  Cross-Sectional  Area 

For  the  dome  the  cross-sectional  area  is 

A  =  R2  cos"1  (•L^)  ~  (  R  *  h4  )  (  2  R  h4  -  h2  )1/2 

For  the  sides  the  area  is  simply 
A  =  2  r1  (  h2  +  h3  ) 

The  total  sross-sectional  area  is  then 

A  =  3.38  x  103  ft3  (  3.14  x  102  m2  ) 
xs 
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APPENDIX  C 


DERIVATION  OF  THE  DISK-SOURCE 
FLUX  EQUATION  (Eg.  4.9)  [27,29] 


For  a  point  source  of  strength  S  photons/sec  shielded  by  a  slab 
b  mean  free  paths  thick  the  flux  reaching  point  P  (see  Fig.  C.l)  at 
a  distance  of  x  cm  is: 


B  S  -b 

K  «  - 2  6 

Y  4ttx 


photons/cm  -sec 


where  B  =  Point  isotropic  build-up  factor 

(a  function  of  shield  material,  b  and  E) 


Therefore  for  a  disk  of  source  strength  photons/cm  -sec.  the  flux 
at  P  from  the  differential  ring  between  r  and  r  +  dr  is : 


B  SA  (2lrrdr)  _bsece„ 


„  2  2  2 
But  P  =  r  +  x 


therefore  pdp  =  rdr 


Also  sec0"  =  — 
x 


therefore 


BSA  dp  -b 
~ 2  p~  X 


Figure  C.l 

Geometry  for  Disk  Source  Flux  Derivation  [31] 
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Substituting 


and  integrating  over  the  disk. 


Y 


bsec0 


b 


/ 


dt 


Since,  by  definition, 


E1(x) 


the  flux  can  be  expressed  as 


Y 


BS  2 

— t—  [E^ (b)  -  E^(bsecB) ]  photons/cm 


-sec 
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APPENDIX  D 


DERIVATION  OF  THE  CYLINDRICAL  VOLUME  SOURCE 


FLUX  EQUATION  (Eg.  4.19)  [27] 


For  the  self-absorbing  cylindrical  volume  source  shown  in  Fig.  D.l 
the  uncollided  flux  at  point  P  is: 


BS  -y  x  -yy 

i _ v  r  e  s  e  ... 

S'  Ait  V  2  dV 

r 


where 


Distance  from  dV  to  the  cylinder  surface 
(the  self-absorption  distance) 


=  Attenuation  path  length  in  the  shield 


=  Distance  from  dV  to  P 


=  Linear  attenuation  coefficient  in  the  source 
=  Linear  attenuation  coefficient  in  the  shield 


For  the  cylindrical  coordinate  system  z,p,f 


dV  =  dzdpd<p 


v  ^  /  dz  /  pdp  /  v-f -T4* - 

0  0  Op  +  s  +  z-  2spcos‘f> 


It  can  be  seen  from  Fig.  D.l  that 


x  =  (PB  -  PC)  sec  ^ 


where 


PD  =  (z2+  p2+  s2  -  2spcosy3  )] 
PB  (p2+  s2  -  2spcosf  )1/2 
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k  1  7T  b2 

G  (  k.p.uR.b,  )  =  )  dn  f  mdm  /  — x - x - x - 

0  0  On  +  m+  p-  2mpcos  *p 


,  2^  2^  2  ,  ..  ,1/2  . 

exp  {-(n+m+p  -  2mpcos ) 


2  2  2  2  2  2  1  19 

T?  m  ~  mpcosy  +  (  m  +  p  -  2mpcosy  -  m  p  sin  'P  ) 

^sK  2  2  + 

m  +  p  -  2mpcos  f 


The  function  G  (  k,p,ysR,b2  )  has  been  numerically  evaluated  and 

listed  for  selected  values  of  k,  p,  ygR,  and  in  [27],  Note  that  for 

no  self-absorption  y  R  =  0.  The  values  of  the  G  function  listed  in 

s 

Table  A. 10  were  obtained  by  graphing  the  function  as  listed  in  [27] 
and  interpolating  to  the  required  values  of  b^  and  b'^. 


